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Abstract Recent years have seen increasing efforts to 
directly measure some aspects of the general relativis- 
tic gravitomagnetic interaction in several astronomical 
scenarios in the solar system. After briefly overview- 
ing the concept of gravitomagnctism from a theoretical 
point of view, we review the performed or proposed at- 
tempts to detect the Lense-Thirring effect affecting the 
orbital motions of natural and artificial bodies in the 
gravitational fields of the Sun, Earth, Mars and Jupiter. 
In particular, we will focus on the evaluation of the im- 
pact of several sources of systematic uncertainties of 
dynamical origin to realistically elucidate the present 
and future perspectives in directly measuring such an 
elusive relativistic effect. 
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1 Introduction 

The analogy between Newton's law of gravitation and 
Coulomb's law of electricity has been largely investi- 
gated since the nineteenth century, focusing on the 
possibility that the motion of masses could produce 
a magnetic-like field of gravitational origin. For in- 
stance, Holzmiiller (1870) and Tisserand (1872, 1890), 
taking into account the modification of the Coulomb 
law for the electrical charges by Weber (1846), pro- 
posed to modify Newton's law in a similar way, intro- 
ducing in the radial component of the force law a term 
depending on the relative velocity of the two attract- 
ing particles, as described by North (1989) and Whit- 
taker (1960). Moreover, Hcaviside (1894) investigated 
the analogy between gravitation and elcctromagnctism; 
in particular, he explained the propagation of energy in 
a gravitational field in terms of an electromagnetic-type 
Poynting vector. 

Actually, today the term "gravitomagnctism" (CM) 
(Thornc 1988; Rindlcr 2001; Mashhoon 2007) com- 
monly indicates the collection of those gravitational 
phenomena regarding orbiting test particles, process- 
ing gyroscopes, moving clocks and atoms and propa- 
gating electromagnetic waves (Dymnikova 1986; Rug- 
giero and Tartaglia 2002; Schafer 2004, 2009) which, 
in the framework of the Einstein's General Theory of 
Relativity (GTR), arise from non-static distributions 
of matter and energy. In the weak-field and slow 
motion approximation, the Einstein field equations of 
GTR, which is a highly non-linear Lorentz-covariant 
tensor theory of gravitation, get linearized, thus look- 
ing like the Maxwcllian equations of electromagntism. 
As a consequence, a "gravitomagnetic" field B g , in- 
duced by the off-diagonal components goi, i = 1,2,3 
of the spacetime metric tensor related to mass-energy 
currents, does arise. Indeed, bringing together Newto- 
nian gravitation and Lorentz invariance in a consistent 
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ficld-thcorctic framework necessarily requires the intro- 
duction of a "magnetic" -type gravitational field of some 
form (Khan and O'Connell 1976; Bedford and Krumm 
1985; Kolbenstvedt 1988). 

In general, GM is used to deal with aspects of GTR 
by means of an electromagnetic analogy. However, it 
is important to point out that even though the lin- 
earization of the Einstein's field equations produces the 
Maxwell-like equations (the so called " linear perturba- 
tion approach" to GM, see e.g. Mashhoon (2007)), of- 
ten written in the literature including time dependent 
terms, they are, in that case, just formal (i.e., a differ- 
ent notation to write linearized Einstein equations) , as 
the 3-vectors E g and B g (the "gravito-electromagnetic 
fields") showing up therein do not have a clear phys- 
ical meaning. A consistent physical analogy involv- 
ing these objects is restricted to stationary phenom- 
ena only (Clark and Tucker 2000; Costa and Herdeiro 
2008, 2010), that is, actually, the case treated here. 
One may check, for instance, that from the geodesies 
equation the corresponding Lorentz force is recovered - 
to first order in v/c - only for stationary fields (Costa 
and Herdeiro 2008; Bini et al. 2008). Moreover, the 
Maxwell-like equations obtained by linearizing GTR 
have limitations 1 , since they are self-consistent at lin- 
ear order only, which is what we are concerned with in 
this paper; in fact, inconsistencies arise when this fact 
is neglected 2 . 

Far from a localized rotating body with angular mo- 
mentum S the gravitomagnctic field can be written 
as (Thornc et al. 1986; Thornc 1988; Mashhoon et al. 
2001a) 



B g {r) 



G 
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(1) 



where G is the Newtonian gravitational constant, c the 
speed of light in vacuum and f a unit vector along r. 
Eq. (1) equals the field of a magnetic dipole with a 
moment fl' g = GS/c and affects, e.g., a test particle 
moving with velocity v with a non-central acceleration 
(Mashhoon et al. 2001b; Lichtenegger et al. 2006) 



^4gm — —2 



x B„, 



(2) 



which is the cause of three of the most famous and em- 
pirically investigated GM features with which we will 
deal here: the Lense and Thirring (1918) effect, the gy- 
roscope precession (Pugh 1959; Schiff 1960a), and the 
gravitomagnetic clock effect (Mashhoon et al. 2001a). 

1 See, e.g., Misncr ct al. (1973) Section 7.1, box 7.1 and Section 
18.3, and Ohanian and RufBni (1994), Chapter 3. 

2 See, e.g., Tartaglia and Ruggiero (2004). 



Corrections of higher order in v/c to orbital motions of 
pointlikc objects have been studied by Capozziello et 
al. (2009); they may become relevant in stronger gravi- 
tational fields like those occurring in astrophysical sce- 
narios (pulsars and black holes). For recent reviews of 
the Lense-Thirring effect in the astrophysical context 
see, e.g., Stella and Possenti (2009); Schafer (2009). 

GM manifests itself in the gravitational tidal forces 
as well. Actually, a gravito-electromagnetic analogy re- 
lying on exact and covariant equations stems from the 
tidal dynamics of both theories (Costa and Herdeiro 
2008). In both theories it is possible to define elec- 
tric and magnetic-type tidal tensors playing analogous 
physical roles. The electric tidal tensors of electromag- 
nctism are gradients of the electric field, and play in 
the worldline deviation for two neighboring test parti- 
cles (with the same ratio charge to mass) the same role 
as the so called "electric part of the Riemann tensor" 
in the geodesic deviation equation. The magnetic-type 
tidal tensors are, in the electromagnetic case, gradients 
of the magnetic field, and their gravitational counter- 
part is the so-called magnetic part of the Riemann ten- 
sor. Physically, the latter manifests itself, for instance, 
in the deviation from geodesic motion of a spinning 
test particle due to a net gravitational force acting on 
it, in analogy with the electromagnetic force exerted on 
a magnetic dipole. In the framework of the tidal tensor 
formalism, the exact gravitational force on a gyroscope 
(Mathisson 1937; Papapetrou 1951; Pirani 1956) is de- 
scribed by an equation formally identical to the electro- 
magnetic force on a magnetic dipole. Moreover, both 
Maxwell and part of the Einstein equations (the time- 
time and time-space projections) may be expressed ex- 
actly as equations for tidal tensors and sources, and 
such equations exhibit a striking analogy. In particular, 
these equations show explicitly that mass currents gen- 
erate the gravitomagnetic tidal tensor just like currents 
of charge involve the magnetic tidal tensor of electro- 
magnetism (Costa and Herdeiro 2008). At this point it 
is important to make a distinction between the "gravit- 
omagnetic field" causing, e.g., the precession of a gyro- 
scope (see Section 1.2), and the gravitomagnetic tidal 
field causing e.g. the net force on a the gyroscope. The 
gravitomagnctic field itself has no physical existence; it 
is a pure coordinate artifact that can be gauged away 
by moving to a freely falling (non-rotating) frame. For 
instance, it is well known that the spin 4- vector of a gy- 
roscope undergoes Fermi- Walker transport, with no real 
torques applied on it; thus the gyroscope "precession" is 
a non-covariant notion attached to a specific coordinate 
system, which is the reference frame of the distant stars, 
i.e., it does not have a local existence, and can only be 
measured by locking to the distant stars by means of a 
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telescope (Polnarev 1986). The gravitomagnetic tidal 
tensor, instead, describes physical forces, which can be 
locally measured (Polnarev 1986). 

There have been several proposals to detect the GM 
tidal forces originating by the proper angular momen- 
tum S of a central body like the Earth. In principle, 
they can be detected by an orbiting tidal force sensor, 
or "gravity gradiometer" , as suggested by Braginsky 
and Polnarev (1980). In the following years such a con- 
cept was further investigated (Mashhoon and Theiss 
1982), also from the point of view of a practical im- 
plementation in terms of a set of orbiting supercon- 
ducting gravity gradiometers (SGG) (Mashhoon et al. 
1989; Paik 1989). Paik (2008) recently reviewed the 
the feasibility of the gravity gradiometer experiment in 
view of the latest advancements in the field. Accord- 
ing to Paik (2008), the GM field of the Earth could be 
successfully detected by using an orbiting SGG. Such a 
mission would benefit from the technologies already de- 
veloped for GP-B (see Section 1.2). The SGG would be 
launched in a superfluid helium dewar, and the helium 
boil-off gas would be used for drag-free control of the 
satellite. GP-B, like gyros or telescopes, may be used 
for attitude control of the spacecraft. 

Other GM effects which have recently received at- 
tention from the phenomenological point of view are 
those caused by the orbital motion of the Earth-Moon 
system around the Sun. They have nothing to do with 
the proper angular momenta S of the Earth or the Sun 
causing the "intrinsic" GM effects like those previously 
mentioned; in the case of translational mass-energy cur- 
rents it is customarily to speak about "extrinsic" GM 
effects. According to Nordtvedt (1988, 2003), the ex- 
trinsic GM interaction has already been observed with a 
relative accuracy of 1 part to 1000 in comprehensive fits 
of the motions of several astronomical and astrophysical 
bodies like satellites, binary pulsars and the Moon. In 
fact, a debate arose about the ability of the Lunar Laser 
Ranging (LLR) technique (Dickey et al. 1994) of detect- 
ing genuine GM effects (Murphy et al. 2007a; Kopcikin 
2007; Murphy et al. 2007b; Soffel et al. 2008; Ciufolini 
2010). Concerning the Lense-Thirring precessions of 
the lunar orbit induced by the Earth's spin, they may 
remain still undetectable also in the foreseeable future 
because of overwhelming systematic uncertainties (Iorio 
2009a), in spite of the expected mm- level improvements 
in Earth- Moon ranging with LLR (Murphy et al. 2008). 

Several Earth-based laboratory experiments aimed 
to test the influence of the "intrinsic" terrestrial GM 
field on classical and quantum objects and electromag- 
netic waves have been proposed so far (Braginsky et 
al. 1977, 1984; Cerdonio et al. 1988; Ljubicic and Lo- 
gan 1992; Camacho and Ahluwalia 2001; Tartaglia and 



Ruggiero 2002; Iorio 2003a; Stedman et al. 2003; Iorio 
2006a), but they have never been implemented because 
of several technological difficulties in meeting the strin- 
gent requirements in terms of sensitivity and/or accu- 
racy. 

As gravitomagnetic effects are, in general, analyzed 
in the framework of the linearized theory of GTR, they 
have also an important connection with Gravitational 
Waves (GWs) (Iorio and Corda 2009, 2010; Baskaran 
and Grishchuk 2004; Corda 2007b). We recall that the 
data analysis of interferometric detectors has nowadays 
been started, and the scientific community awaits a first 
direct detection of GWs in the near future. For the cur- 
rent status of GW interferometers we refer to Giazotto 
(2008). Thus, the indirect evidence for the existence of 
GWs by Hulse and Taylor (1975) might soon be con- 
firmed. Detectors for GWs will be important for a bet- 
ter knowledge of the Universe (Iorio and Corda 2010) 
and also because the interferometric GWs detection will 
be the ultimate test for GTR or, alternatively, a strong 
endorsement for Extended Theories of Gravity (Corda 
2009). In fact, if advanced projects on the detection 
of GWs improve their sensitivity, allowing the scien- 
tific community to perform a GW astronomy, accurate 
angle- and frequency-dependent response functions of 
interferometers for GWs arising from various theories 
of gravity will permit to discriminate among GTR and 
extended theories of gravity. This ultimate test will 
work because standard GTR admits only two polariza- 
tions for GWs, while in all extended theories there exist 
at least three polarizations states; see Corda (2009) for 
details. On the other hand, the discovery of GW emis- 
sion by the compact binary system composed of two 
Neutron Stars PSR1913+16 (Hulse and Taylor 1975) 
represents, for scientists working in this research field, 
the definitive thrust allowing to reach the extremely so- 
phisticated technology needed for investigating in this 
field of research (Iorio and Corda 2010). GWs are a 
consequence of Einstein's GTR (Einstein 1915), which 
presupposes GWs to be ripples in the spacetime curva- 
ture traveling at light speed (Einstein 1916, 1918; Iorio 
and Corda 2010). The importance of gravitomagnetic 
effects in the field of a GW has been emphasized by 
Baskaran and Grishchuk (2004). For a complete review 
of such a topic, see Iorio and Corda (2010). Recently, 
the analysis has been extended to gravitomagnetic ef- 
fects in the field of GWs arising by Scalar Tensor Grav- 
ity too (Corda et al. 2010; Iorio and Corda 2010). 

1.1 The Lense-Thirring effect 

After the birth of the Einstein's Special Theory of Rela- 
tivity (STR) in 1905, the problem of a "magnetic" -type 
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component of the gravitational field of non-static mass 
distributions was tackled in the framework of the search 
for a consistent relativistic theory of gravitation (Ein- 
stein 1913). 

With a preliminary and still incorrect version of 
GTR, Einstein and Besso in 1913 (Klein et al. 1955) 
calculated the node precession of planets in the field 
of the rotating Sun; the figures they obtained for Mer- 
cury and Venus were incorrect also because they used 
a wrong value for the solar mass. Soon after GTR was 
put forth by Einstein in 1915, de Sitter (1916) worked 
out the corresponding shifts of the planet's perihelia 
for ecliptic orbits due to the rotation of the Sun; how- 
ever, his result for Mercury (—0.01 arcseconds per cen- 
tury) was too large by one order of magnitude because 
he assumed a homogenous and uniformly rotating Sun. 
In 1918 Thirring (1918a) analyzed in a short article 
the formal analogies between the Maxwell equations 
and the linearized Einstein equations. Later, Thirring 
(1918a, b, 1921) computed the centrifugal and Coriolis- 
like gravitomagnetic forces occurring inside a rotating 
massive shell. Lense and Thirring (1918) 3 worked out 
the gravitomagnetic effects on the orbital motions of 
test particles outside a slowly rotating mass; in partic- 
ular, they computed the gravitomagnetic rates for the 
two satellites of Mars (Phobos and Deimos), and for 
some of the moons of the giant gaseous planets. They 
found for the longitude of the ascending node 57 and the 
argument of pericenter uj a pro- and retrograde preces- 
sion, respectively, according to 



fl hT = 



2GS 



c 2 a 3 (l - e 2 ) 3 / : 



6GScosI 

c 2 a 3( 1 _ e 2)3/2' 



(3) 



where a, e, / are the semimajor axis, the eccentricity 
and the inclination of the test particle's orbital plane 
to the central body's equator, respectively Later, 
the Lense-Thirring effect was re-derived in a number 
of different approaches; sec, e.g., Bogorodskii (1959); 
Zel'dovich and Novikov (1971); Barker and O'Conncll 
(1974); Landau and Lifshitz (1975); Ashby and Alli- 
son (1993); Iorio (2001a); Lammerzahl and Neugebauer 
(2001); Chashchina et al. (2009). In the following we 
give some more physical insights about the derivation 
and the characteristics of cq. (l)-cq. (2) and, conse- 

3 However, in August 1917 Einstein (Schulmann et al. 1998) wrote 
to Thirring that he calculated the Coriolis-type field of the ro- 
tating Earth and Sun, and its influence on the orbital elements 
of planets (and moons). A detailed history of the formulation of 
the so-called Lense-Thirring effect has recently been outlined by 
Pfister (2007); according to him, it would be more fair to speak 
about an Einstein- Thirring-Lcnse effect. 

4 In fact, Lense and Thirring (1918) considered the longitude of 
pericenter to = $1 + u>. 



quently, of eq. (3). In the standard parameterized post- 
Newtonian 5 (PPN) framework for an isolated, weakly 
gravitating and slowly rotating body, the PPN space- 
time metric coefficients involving its angular momen- 
tum are the off-diagonal terms 



9oi = 
with 

V 



(7+1) 



GS xf 



V,, 1 = 1,2,3, 



(4) 



(5) 



In eq. (4) 7 is the PPN parameter accounting for the 
spatial curvature caused by a unit mass; in GTR 7 = 1. 
Thus, in the PPN approximation the resulting Christof- 
fel symbols including Vj, i — 1,2,3 and entering the 
geodesic equation of motion are (Soffel 1989) 



1 r 



0j 



(7 + 1) f dVj _ dVj_ 
c 2 V dxi dx l 



i,j = 1,2,3; 



(6) 



they give rise to eq. (l)-eq. (2) for 7=1. It is in- 
teresting and useful to emphasize the complete analogy 
of eq. (2) with the Lorentz force of electromagnctism 
provided that the factor —2 entering eq. (2) can be 
formally thought as a gravitomagnetic chargc-to-mass 
ratio q g /m. As originally done by Lense and Thirring 
(1918), the orbital precessions of the node f2 and the 
pericenter u> of eq. (3) can be straightforwardly worked 
out with, e.g., the Gauss equations for the variation of 
the Keplerian orbital elements (Soffel 1989; Roy 2005) 
by treating the gravitomagnetic force eq. (2) as a small 
perturbation of the Newtonian monopole. The Gauss 
equations for ft and cj are (Soffel 1989; Roy 2005) 



~dt 



1 



ssiniyl — e 2 



-.A 



v sinu, 



(7) 



dt nae 



-A r cos / + A T (1 + £J sin/ 



cos/f , 



(8) 



5 The post-Newtonian formalism allows to approximate the non- 
linear Einstein field equations for weak fields and slow motions 
in terms of the lowest order deviations from the Newton's theory. 
In the framework of the parameterized post-Newtonian formal- 
ism (Eddington 1922; Ni 1972; Nordtvcdt 1968, 1969; Will 1971), 
such departures from classical gravity are expressed in terms of 
a set of ten parameters to discriminate between competing met- 
ric theories of relativistic gravity. The modern notation is due 
to Will and Nordtvcdt (1972). For a comprehensive review, see 
Will (1993). 
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where n = y/ GM/a 3 is the Keplerian mean motion, 
u = u; + / is the argument of latitude in which / is 
the true anomaly, and A r , A T , A v are the projections of 
the perturbing acceleration onto the radial, transverse 
and normal components of an orthonormal frame co- 
moving with the test particle. In the case of eq. (2) 
these projections correspond to (Soffel 1989) 



<LT 



A^ T = — 7?e cos / sin /, 



r\ cos 1(1 + ecos/), 



A,/ = r) sin 1(1 + ecos/) 



(9) 



(10) 



2 sin u + e 



sin / cos u 
1 + e cos / 



with 



Xn 



a 2 (l - e 2 ) 7 /: 



■(1 + ecos/) 3 , 



and 



A 



2GS 



(11) 



(12) 



(13) 



It is immediately seen from eq. (9)-eq. (11) that for 
polar orbits (i.e. I = 90 deg), both the radial and the 
transverse components vanish, contrary to the out-of- 
plane component. This fact can easily be inferred also 
from the dipolar field structur and the Lorentz-type 
force, i.e. eq. (l)-eq. (2). Indeed, if S lies in the or- 
bital plane, the latter one contains B g as well, so that, 
according to eq. (2), Aqm is entirely out-of-plane. In 
that case there is also no pcricenter precession, as in- 
dicated by eq. (3) since all terms in eq. (8) vanish for 
I = 90 deg. By analogous reasoning it is easy to fig- 
ure out that for equatorial orbits (i.e. I — deg) Aqm 
completely lies in the orbital plane which is in agree- 
ment with eq. (9)-eq. (11). If the orbit is circular in 
addition, Aqm, which is perpendicular to S and v, be- 
comes entirely radial, as confirmed by eq. (9)-eq. (10). 
Finally, since the factor sin J cancels in eq. (7) and eq. 
(11) the Lense-Thirring node precession is independent 
of /, as stated by eq. (3). 

1.2 The gyroscope precession 

Another well known GM effect consists of the preces- 
sion of a gyroscope moving in the field of a slowly 
rotating body. It was worked out in 1959 by Pugh 
(1959) and in 1960 by Schiff (1960a,b,c) on the basis 
of the Mathisson (1937) and Papapetrou (1951) equa- 
tion and became known as the Schiff effect. More recent 
derivations, based on a quantum mechanical approach 



to gravitation, can be found in Barker and O'Connell 
(1970, 1972). 

In order to yield a direct, physical insight of such 
a phenomenon, let us recall the following basic facts 
of Maxwellian electromagnetism. A charged spinning 
particle with electric charge q, mass m and spin a has 
a magnetic moment 
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2mc 



(14) 



so that it precesses in an external magnetic field B ac- 
cording to 



da 



ft x B 



I 1 



\2mc 

with the precessional frequency 



B x a = D x a 



O = - 



<1 



2mc 



B. 



(15) 



(16) 



Moving now to the weak-field and slow-motion ap- 
proximation of linearized gravitomagnetism and recall- 
ing that the gravitomagnetic charge-to-mass ratio is 
q g /m = —2, according to eq. (14) a spinning parti- 
cle like a gyroscope is endowed with a gravitomagnetic 
dipole moment 



Ms 



a 

* 

c 



(17) 



so that it undergoes a precession during its motion in 
an external gravitomagnetic field B g with frequency 



6 =^ 

l-t g — 

C 



(18) 



Since the GM-ficld of a distant rotating astronomical 
body is given by eq. (1), the precession frequency be- 
comes 



G 



S -3 S-f)f 



(19) 



If the gyroscope moves along an equatorial orbit 6 , then 
da G 



dt 



-S xa 



(20) 



and there is no precession of the orbit in case the two 
spins arc aligned. If the orbit of the gyroscope is polar, 
then, by choosing a plane {xz} reference frame with, 
say, 



S = S 



(21) 



3 For simplicity, we will assume it circular. 
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and 

f = cos nt x + sin nt i 



(22) 



it is easy to show that the averaged precession vanishes 
if a and S are parallel or antiparallcl, i.e. for 



a = ±er z. 



(23) 



Finally, if a is orthogonal to S and to the orbital plane, 
i.e. 



a = ±<jy, (24) 
then there is a net spin precession given by 
ds\ G 



dt / 2c 2 r 



■S x a. 



(25) 



Soon after the formulation of the Schiff effect, in 
1961 Fairbank and Schiff (1961) submitted to NASA 
a proposal for a dedicated space-based project aimed 
to directly measure the precession of eq. (25) in a ded- 
icated, controlled experiment. Such an extraordinary 
and extremely sophisticated mission, later named Grav- 
ity Probe B (GP-B) (Evcritt 1974; Evcritt ct al. 2001), 
consisted of a drag-free, liquid helium-cooled spacecraft 
moving in a polar, low 7 orbit around the Earth and car- 
rying onboard four superconducting gyroscopes whose 
GM precessions of 39 milliarcseconds per year (mas 
yr _1 in the following) should have been detected by Su- 
perconducting Quantum Interference Devices (SQUID) 
with an expected accuracy of 1% or better. It took 
43 years to be implemented before GP-B was finally 
launched on 20 April 2004; the science data collec- 
tion lasted from 27 August 2004 to 29 September 2005, 
while the data analysis is still ongoing (Conklin et al. 
2008; Everitt et al. 2009). It seems that the final accu- 
racy obtainable will be less than initially expected be- 
cause of the occurrence of unexpected systematic errors 
(Muhlfelder et al. 2009; Keiser et al. 2009; Silbergleit et 
al. 2009). At present, the GP-B team reports 8 evidence 
of the gravitomagnetic spin precessions as predicted by 
GTR, with a statistical error of « 14% and systematic 
uncertainty of w 10%. 

In 1975 Haas and Ross (1975) proposed to measure 
the angular momenta of the Sun and Jupiter by exploit- 
ing the Schiff effect with dedicated spacecraft-based 
missions, but such a proposal was not carried out so 
far. 

7 Its altitude was 642 km. 

8 See on the WEB: http://einstein.stanford.edu/ 



1.3 The gravitomagnetic clock effect 

Zel'dovich (1965) discovered that the gravitational field 
of a rotating mass of radius R and angular momentum 
S splits the line emitted by an atom with frequency fo 
into two components with opposite circular polariza- 
tions and frequencies fo ± Af, where 



Af: 



2GS 
c 2 R 3 



(26) 



for an electromagnetic signal emitted on the body's sur- 
face at the pole. It is a gravitational analog of the 
electromagnetic Zeeman effect and, as such, it does not 
depend on the specific properties of the system emit- 
ting the electromagnetic radiation, thus being the same 
for an atom or a molecule and independent from the 
emitted frequency. A similar effect for the orbital mo- 
tion of test particles around a spinning body of mass 
M and angular momentum S was discovered later by 
Vladimirov et al. (1987). They noted that in the equa- 
torial plane of a rotating mass the gravitomagnetic force 
becomes purely radial (if the motion occurs at a con- 
stant distance r), and acts as a supplement to the cen- 
tripetal Newtonian monopole; it is clearly elucidated 
by eq. (9)-eq. (10). Such an additional term, may be 
either positive, i.e. directed outward so that it weakens 
the overall gravitational force, if the particle's motion 
coincides with the direction of rotation of the central 
body, or negative, thus enhancing the force of gravity, 
if the motion is opposite to the rotation. As a conse- 
quence, the orbital period of a particle moving along a 
circular and equatorial orbit in the Lensc-Thirring met- 
ric is longer in the first case and shorter in the second 9 . 
This fact contradicts the idea of frame-dragging where 
it is conceived that a moving object is "dragged" by 
spacetime which in turn is "twisted" by the rotation of 
the central mass. If this were the case, the "dragged" 
test particle should move faster when co-rotating with 
the central body and should thus have a shorter period. 
Interestingly, such a change in the period of the parti- 
cle's orbit is a universal and structurally very simple 
quantity, given by 



GM 



±2tt 



S 
Mc 2 ' 



(27) 



It is amazing that eq. (27) is independent of both the 
Newtonian gravitational constant G and the orbital ra- 
dius. Cohen and Mashhoon (1993) suggested to con- 
sider the difference between the orbital periods of two 

9 Rcmcmber that the period of a pendulum's oscillations decreases 
as the force acting on it increases 
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counter-orbiting clocks moving around a rotating as- 
tronomical body along circular and equatorial orbits 
because it cancels the common Keplerian terms and 
enhances the gravitomagnetic ones by adding them up. 
This is the so-called gravitomagnetic clock effect (Mash- 
hoon ct al. 1999, 2001a). Theiss (1985) worked out the 
case of circular orbits with arbitrary inclinations show- 
ing that the time difference decreases with increasing 
inclination; for a polar orbit the effect vanishes as it is 
expected because of symmetry. The general case for ar- 
bitrary values of the eccentricity and the inclination was 
treated by Mashhoon et al. (2001b) and Lichtcncggcr 
and Iorio (2007). If one considers two satellites orbiting 
a slowly rotating mass M in opposite directions, their 
common initial position in the orbital plane is given by 
the argument of latitude u$ = ujq + /q. In taking the 
difference of the sidereal periods of the satellites, the 
gravitoelectric perturbations cancel, leaving 



AT 



sid 



47r5'cos/ 
c 2 M 



2 (2 — tan 2 J cos 2 u ) 



(1 + ecos/o) 2 



(28) 



which reduces to the difference of the GM corrections 
Tqm for counter-orbiting particles of (27) in case of e = 
1 = 0. It is interesting to note that this clock effect can 
reveal a relatively large value by a careful choice of the 
initial parameters (Lichtenegger and Iorio 2007). 

Based on various approaches (e.g., analogies with 
electromagnetism, spacetimc geometric properties), 
derivations of the gravitomagnetic clock effect for the 
circular and equatorial cases can be found in (You 1998; 
Iorio ct al. 2002b; Tartaglia 2000a). Gronwald ct al. 
(1997) proposed to detect the gravitomagnetic clock 
effect with a space-based mission - dubbed Gravity 
Probe C - in the gravitational field of the Earth where 
such an effect would be as large as 10~ 7 s. In view 
of the challenging difficulties of implementing such a 
demanding experiment, a number of studies were per- 
formed to mainly investigate the impact of several com- 
peting dynamical effects acting as sources of insidious 
systematical uncertainty (Lichtenegger et al. 2000; Io- 
rio 2001b, c; Iorio and Lichtenegger 2005; Lichtenegger 
et al. 2006). Tartaglia (2000b) preliminarily looked at 
the possibilities offered by other Solar System scenarios 
finding that, in principle, the less unfavorable situation 
occurs for the Sun and Jupiter. 



1.4 Motivations for attempting to directly measure 
the Lense-Thirring effect 

GTR is a basic pillar of our knowledge of Nature 
since it currently represents our best theory of gravi- 



tation, which is one of the four fundamental interac- 
tions governing the physical world. The simplicity, the 
internal coherence and the mathematical elegance of 
GTR are remarkable, but the level of its empirical cor- 
roboration, although certainly satisfactorily up to now 
(Will 2006), is not comparable to that of the other theo- 
ries describing the remaining fundamental interactions. 
This applies both to the number of successfully tested 
predictions of GTR as well as to the level of accuracy 
reached. Notably the level of empirical corroboration 
of GM, which is a constitutive, fundamental aspect of 
GTR is to date extraordinarily poor. It is therefore 
desirable to expand and strengthen the empirical ba- 
sis of the theory by testing as many diverse aspects 
and predictions as possible using different methods and 
techniques. It is also necessary, however, to devise ob- 
servational/experimental tests in those extreme regimes 
in which the theory in its currently accepted form is be- 
lieved to experience failures. This is particularly impor- 
tant in order to gain possible hints on a future quantum 
gravity theory which should combine both gravitational 
and quantum phenomena. 

Concerning the GP-B mission, aimed to test a well de- 
fined GM prediction using the gravitational field of the 
Earth, although different teams may have the possi- 
bility to repeat the analysis of the currently available 
data record with various approaches and techniques, 
the results of the mission are likely doomed to remain 
unique since it will be impossible to replicate the en- 
tire experiment in any foreseeable future. For the mo- 
ment, its level of accuracy is more or less comparable 
to or better than that reached in the non-dedicated 10 
tests of the Lense-Thirring effect with the terrestrial 
LAGEOS satellites. Thus, it becomes of the utmost im- 
portance not only to reliably assess the total accuracy 
in such attempts but also to look for other possibilities 
offered by different astronomical scenarios by exploiting 
future planned/proposed spacecraft-based missions and 
expected improvements in ranging techniques. Such an 
effort has the merit, among other things, to realistically 
establish the limits which may likely be reached in such 
an endeavor. Moreover, as a non-negligible by-product, 
the knowledge of several classical effects, regarded in 
the present context as sources of unwanted systematic 
biases but interesting if considered from different points 
of view, will turn out to be greatly improved by the 
tireless efforts towards the measurement of such a tiny 
relativistic feature of motion. Last but not least, it 
should be recalled that reaching a satisfying level of 
knowledge of GM has important consequences in the 

10 LAGEOS and LAGEOS II were originally launched for other 
purposes. 
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study of extreme astrophysical scenarios in which, as 
we know, GM may play a very important role (Thorne 
et al. 1986; Thorne 1988; Williams 1995, 2004; Arvani- 
taki and Dubovsky 2010). 

Another source of motivation to exploit gravitomag- 
netic effects is their possible relation with Mach's prin- 
ciple. While GTR as a whole - despite its name - 
appears not to fulfill Machian expectations of a de- 
scription of motion with only relative concepts, some 
special GM phenomena seem to be in accordance with 
Machian ideas. In particular dragging effects are con- 
sidered to be the most direct manifestations of Mach's 
principle in general relativity because they indicate that 
local incrtial frames are at least partially determined by 
the distribution and currents of mass-energy in the uni- 
verse. Thus the study of gravitomagnctism may provide 
a deeper insight into the presumed intimate connection 
between inertial properties and matter. 
2 Measuring the Lense-Thirring effect in the 
gravitational field of the Earth 

Soon after the dawn of the space age with the launch 
of Sputnik in 1957 Soviet scientists proposed to di- 
rectly test the general relativistic Lense-Thirring effect 
with artificial satellites orbiting the Earth. In par- 
ticular, Ginzburg (1957, 1959, 1962) proposed to use 
the perigee of a terrestrial spacecraft in a highly el- 
liptic orbit, while Bogorodskii (1959) considered also 
the node. Yilmaz (1959), aware of the aliasing effect 
of the much larger classical precessions induced by the 
non-sphericity of the Earth, proposed to launch a satel- 
lite in a polar orbit to cancel them. About twenty 
years later, Van Patten and Everitt (1976a); van Patten 
and Everitt (1976b) suggested to use a pair of drag- 
free, counter-orbiting terrestrial spacecraft in nearly 
polar orbits to detect their combined Lense-Thirring 
node precessions. Almost contemporaneously, Cugusi 
and Proverbio (1977, 1978) suggested to use the pas- 
sive geodetic satellite LAGEOS, in orbit around the 
Earth since 1976 and tracked with the Satellite Laser 
Ranging (SLR) technique (Degnan 1985), along with 
the other existing laser-ranged targets to measure the 
Lense-Thirring node precession. About ten years later, 
Ciufolini (1986) proposed a somewhat simpler version 
of the van Patten-Everitt mission consisting of looking 
at the sum of the nodes of LAGEOS and of another 
SLR satellite to be launched in the same orbit, apart 
from the inclination which should be switched by 180 
deg in order to minimize the competing classical pre- 
cessions due to the centrifugal oblateness of the Earth. 
Iorio (2003b) showed that such an orbital configura- 
tion would allow, in principle, to use the difference of 
the perigees as well. Test calculations were performed 



by Ciufolini et al. (1996, 1997a) with the LAGEOS 
and LAGEOS II satellites, according to a strategy by 
Ciufolini (1996) involving the use of a suitable linear 
combination of the nodes 57 of both satellites and the 
perigee u> of LAGEOS II in order to remove the impact 
of the first two multipoles of the non-spherical gravi- 
tational potential of the Earth. Latest tests have been 
reported by Ciufolini and Pavlis (2004), Ciufolini et al. 
(2006), Ciufolini (2007), Ciufolini et al. (2009), Lucch- 
esi (2007a), and Ries et al. (2008); Ries (2009) with 
only the nodes of both satellites according to a combi- 
nation of them explicitly proposed by Iorio and Morca 
(2004). The total uncertainty reached is still a matter 
of debate (Iorio 2005a; Ciufolini and Pavlis 2005; Luc- 
chesi 2005; Iorio 2006b, 2007d, 2009c, 2010a; Ciufolini 
et al. 2009) because of the lingering uncertainties in the 
Earth's multipoles and in how to evaluate their biasing 
impact; it may be as large as sa 20 — 30% according 
to conservative evaluations (Iorio 2005a, 2006b, 2007d, 
2009c, 2010a), while more optimistic views (Ciufolini 
and Pavlis 2004, 2005; Ciufolini et al. 2006; Ries et al. 
2008; Ries 2009; Ciufolini et al. 2009) point towards 
ps 10 - 15%. 

2.1 The LAGEOS-LAGEOS II tests 

LAGEOS (Cohen and Smith 1985) was put into or- 
bit in March 1976, followed by its twin LAGEOS II 
(Zerbini 1989) in October 1992; they are passive, spher- 
ical spacecraft entirely covered by retroreflectors which 
allow for their accurate tracking through laser pulses 
sent from Earth-based ground stations. They orbit 
at altitudes of about 6000 km in nearly circular paths 
markedly inclined to the Earth's equator; see Table 1 
for their orbital geometries and Lense-Thirring node 
precessions 11 . The corresponding linear shifts amount 
to about 1.7 m yr" 1 in the cross-track direction 12 at 
the LAGEOS altitudes. 

Since earlier studies (Bogorodskii 1959; Cugusi and 
Proverbio 1978), researchers were aware that a major 
source of systematic errors is represented by the even 



n Cugusi and Proverbio (1978) quoted 40 mas yr -1 for the Lense- 
Thirring node precession of LAGEOS by modeling the Earth as 
a spinning homogeneous sphere. The correct value for the Lense- 
Thirring node precession of LAGEOS was obtained by Ciufolini 
(1986). 

12 A perturbing acceleration like j4qm i s customarily projected 
onto the radial f, transverse f and cross-track directions of 
an orthogonal frame comoving with the satellite (Soffcl 1989); 
it turns out that the Lense-Thirring node precession affects 
the cross-track component of the orbit according to Aj^lt ~ 
asin/Afi LT (cq. (A65), p. 6233 in (Christodoulidis et al. 1988)). 
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Table 1 Orbital parameters and Lense-Thirring node pre- 
cessions of LAGEOS, LAGEOS II, LARES and GRACE for 
S @ = 5.86 x 10 33 kg m 2 s" 1 (McCarthy and Petit 2004). 
The semimajor axis a is in km, the inclination I is in deg, 
and the Lense-Thirring rate J7lt is in mas yr _1 . 



Satellite 


a 


e 


I 




LAGEOS 


12270 


0.0045 


109.9 


30.7 


LAGEOS II 


12163 


0.014 


52.65 


31.5 


LARES 


7828 


0.0 


71.5 


118.1 


GRACE 


6835 


0.001 


89.02 


177.4 



(I = 2, 4, 6, ... ) zonal (m = 0) harmonic coefficients 13 
Je , I = 2, 4, 6 of the multipolar expansion of the classical 
part of the terrestrial gravitational potential, account- 
ing for its departures from spherical symmetry due to 
the Earth's diurnal rotation, induce competing secular 
precessions 14 of the node and the perigee of satellites 
(Kaula 1966) whose nominal sizes are several orders of 
magnitude larger than the Lense-Thirring ones. They 
cannot be removed from the time series of data without 
affecting the Lense-Thirring pattern itself as well. The 
only thing that can be done is to model such a corrupt- 
ing effect as most accurately as possible and assessing 
the impact of the residual mismodclling on the mea- 
surement of the gravitomagnetic effect. In the case of 
the node, the secular precessions induced by the even 
zonals of the geopotential can be written as 



1=2 



(29) 



where the coefficients (l.£,£ = 2,4,6, .. . depend on the 
parameters of the Earth {GM and the equatorial ra- 
dius R) and on the semimajor axis a, the inclination 
I and the eccentricity e of the satellite. For example, 
for £ = 2 the largest precession is due to the first even 
zonal harmonic J2 which, in the small eccentricity ap- 
proximation valid for geodetic satellites, is given by 



■ 3 {R®\ cos/J 2 

^ = -2 n { — ) (1^2)2 ' 



(30) 



The mismodelling in the geopotential-induced preces- 
sions can be written as 



JQgeopot < \ft„ 



1=2 



SJi, 



(31) 



13 The relation among the even zonals Ji and the normalized 
Stokes gravity coefficients Ciq, which arc customarily determined 
in the global Earth's gravity solutions, is Ji = —\J21 + 1 CfQ. 

14 Also the mean anomaly M experiences secular precession due 
to the even zonals, but it is not involved in the measurement of 
the Lense-Thirring effect with the LAGEOS satellites. 



where 5 Ji represents our uncertainty in the knowledge 
of the even zonals Jg . The coefficients VLj of the aliasing 
classical node precessions (Kaula 1966) J7s e °P ot induced 
by the even zonals have been analytically worked out 
up to £ = 20 in the small eccentricity approximation 
by, e.g. Iorio (2003c). 

The three-elements combination used by Ciufolini et 
al. (1996) allowed for removing the uncertainties in J 2 
and J 4 . In Ciufolini et al. (1998a) a « 20% test was 
reported by using the 15 EGM96 Earth gravity model 
(Lemoine et al. 1998); subsequent detailed analyses 
showed that such an evaluation of the total error bud- 
get was overly optimistic in view of the likely unre- 
liable computation of the total bias due to the even 
zonals (Iorio 2003c; Ries et al. 2003a, b). An analo- 
gous, huge underestimation turned out to hold also for 
the effect of the non-gravitational perturbations (Mi- 
lani et al. 1987) like the direct solar radiation pres- 
sure, the Earth's albedo, various subtle thermal effects 
depending on the physical properties of the satellites' 
surfaces and their rotational state (Inversi and Vespe 
1994; Vespe 1999; Lucchesi 2001, 2002, 2003, 2004; Luc- 
chesi et al. 2004; Ries et al. 2003a), which the perigees 
of LAGEOS-like satellites are particularly sensitive to. 
As a result, the realistic total error budget in the test 
reported in (Ciufolini et al. 1998a) might be as large as 
60 — 90% or (by considering EGM96 only) even more. 

The observable used by Ciufolini and Pavlis (2004) 
with the EIGEN-GRACE02S model (Reigber et al. 
2005) and by Ries et al. (2008) with other more re- 
cent Earth gravity models was a linear combination 16 
of the nodes of LAGEOS and LAGEOS II 



LAGEOS II 



(32) 



which was explicitly computed by Iorio and Morea 
(2004) following the approach put forth by Ciufolini 
(1996). Here, 



Cl 



j^LAGEOS 
^LAGEOS II ' 



(33) 



and by means of eq. (30) we find 



Cl 



COS /LAGEOS ( 1 — e LA GEOS II \ f ulageos II 



aLAGEOS 



15 Contrary to the subsequent models based on the dedicated 
satellites CHAMP (http://www-app2.gfz-potsdam.de/pbl/op/ 
champ/) and GRACE (http://www-app2.gfz-potsdam.de/pbl/ 
op/gracc/index_GRACE.html), EGM96 relics upon multidccadal 
tracking of SLR data of a constellation of geodetic satellites in- 
cluding LAGEOS and LAGEOS II as well; thus the possibility of 
a sort of a — priori 'imprinting' of the Lense-Thirring effect itself, 
not solved-for in EGM96, cannot be neglected. 

16 See also (Pavlis 2002; Ries et al. 2003a,b). 
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(34) 

where the values of Table 1 approximately yield 
ci = 0.544. The Lense-Thirring signature of eq. (32) 
amounts to 47.8 mas yr _1 . The combination of eq. (32) 
allows, by construction, to remove the aliasing effects 
due to the static and time-varying parts of the first 
even zonal J2. The nominal (i.e. computed with the 
estimated values of Ji, i = 4,6,...) bias due to the 
remaining higher degree even zonals would amount to 
about 10 5 mas yr _1 ; the need for a careful and reliable 
modeling of such an important source of systematic bias 
is thus quite apparent. Conversely, the nodes of the 
LAGEOS-type spacecraft are directly affected by the 
non-gravitational accelerations at a « 1% level of the 
Lense-Thirring effect (Lucchesi 2001, 2002, 2003, 2004; 
Lucchesi et al. 2004). For a comprehensive, up-to-date 
overview of the numerous and subtle issues concerning 
the measurement of the Lense-Thirring effect see, e.g., 
(Iorio 2007a). 

2.1.1 Conservative evaluation of the impact of the 
mismodelling in the even zonal harmonics 

A common feature of all the competing evaluations so 
far published is that the systematic bias due to the 
static component of the geopotential was always cal- 
culated by using the released (more or less accurately 
calibrated) sigmas aj e of one Earth gravity model so- 
lution at a time for the uncertainties 8J1. Thus, it was 
said that the model X yields a x% error, the model Y 
yields a y% error, and so on. 

Since a trustable calibration of the formal, statistical 
uncertainties in the estimated zonals of the covariance 
matrix of a global solution is always a difficult task to 
be implemented in a reliable way, a much more realistic 
and conservative approach consists, instead, of taking 
the difference 17 



AJ i = \J i (X)-3 e (Y)\, £ = 2,4,6,.. 



(35) 



of the estimated even zonals for different pairs of Earth 
gravity field solutions as representative of the real un- 
certainty SJe in the zonals (Lerch et al. 1994). In 
Tables 2-13 we present our results for the most re- 
cent GRACE-based models released so far by differ- 
ent institutions and retrievable in the Internet at 18 
http://icgem.gfz-potsdam.de/ICGEM/. The models 
used are EIGEN-GRACE02S (Reigber et al. 2005) 

17 See Fig. 5 of Lucchesi (2007b) for a comparison of the estimated 
C40 in different models. 

18 L. I. thanks M Watkins (JPL) for having provided me with the 
even zonals and their sigmas of the JEM01-RL03B model. 



from GFZ (Potsdam, Germany), GGM02S (Tapley 
ct al. 2005) and GGM03S (Tapley et al. 2007) from 
CSR (Austin, Texas), ITG-Grace02s (Mayer-Giirr et 
al. 2006), ITG-Grace03s (Mayer-Giirr 2007) and ITG- 
Grace2010s (Mayer-Giirr et al. 2010) from IGG (Bonn, 
Germany), JEM01-RL03B from JPL (NASA, USA), 
AIUB-GRACE01S (Jaggi et al. 2008) and AIUB- 
GRACE02S (Jaggi et al. 2009) from AIUB (Switzer- 
land). This approach was explicitly followed also, e.g., 
by Milani et al. (1987) with the GEM-L2 and GEM 9 
models, and by Ciufolini (1996) with the JGM3 and 
GEMT-2 models. Note that we do not consider mod- 
els including data from CHAMP, LAGEOS itself 19 and 
Earth-based data. In Tables 2-13 we quote both the 
sum X^=4 ft °f the absolute values of the individual 
mismodcllcd terms (denoted by SAV) 



ft 



^LAGEOS 



C\Ct 



LAGEOS II 



(36) 



and the square root of the sum of their squares 

J J2e=4 fe (RSS); in both cases we normalized them 
to the combined Lense-Thirring total precession of 47.8 
mas yr _1 . 

The systematic bias evaluated with a more realistic 
approach is about 3 to 4 times larger than the one ob- 
tained by using only a single particular model. The 
scatter is still quite large and far from the 5 — 10% 
claimed in Ciufolini and Pavlis (2004). In particular, 
it appears that J4, Jq, and to a lesser extent Js, which 
are just the most relevant zonals for us because of their 
impact on the combination of eq. (32) , are the most un- 
certain ones, with discrepancies AJi between different 
models, in general, larger than the sum of their sigmas 
Gj t , whether calibrated or not. 

Such an approach has been criticized by Ciufolini 
et al. (2009) by writing that one should not compare 
models with different intrinsic accuracies. Moreover, 
Ciufolini et al. (2009) claim that our method would be 
exactly equivalent to compare a modern value of the 
Newtonian gravitational constant G, accurate to 10 -5 , 
to the earlier results obtained in the 18 th century, accu- 
rate to 10 -2 , and conclude that the present-day accu- 
racy would be wrong by a factor 1000. Such criticisms 
are incorrect for the following reasons. According to 
CODATA 20 , the present-day relative accuracy in G is 
1.0 x 10~ 4 = 0.01%, not 0.0015%, corresponding to 



19 It is just one of the devices with which the Lense-Thirring ef- 
fect is measured: using Earth's gravity models including its data 
would yield a-priori "imprint" of GTR itself. Note that the lat- 
est models by GFZ are unsuitable for our purposes because they 
make use of LAGEOS data. 

20 See http://physics.nist.gov/cuu/Constants/ on the WEB. 
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Table 2 Impact of the mismodelling in the even zonal harmonics on f t = p^ AGEOS + cif2 L £ AGEOS 11 AJ e , i = 4, . . . , 20, 
in mas yr -1 . Recall that Jg = — \/2£ + 1 Ciq; for the uncertainty in the even zonals we have taken here the difference 
ACeo = \ Ca > - C^cPI between the model X=GGM02S (Tapley et al. 2005) and the model Y=ITG-Grace02s (Mayer-Gurr 
et al. 2006). GGM02S is based on 363 days of GRACE-only data (GPS and intersatellite tracking, neither constraints nor 
regularization applied) spread between April 4, 2002 and Dec 31, 2003. The a are formal for both models. ACm are always 
larger than the linearly added sigmas, apart from £ = 12 and £ = 18. Values of ft smaller than 0.1 mas yr -1 have not been 
quoted. The Lense-Thirring precession of the combination of eq. (32) amounts to 47.8 mas yr -1 . The percent bias 5/j, have 
been computed by normalizing the linear sum of fi, £ = 4, . . . , 20 (SAV) and the square root of the sum of ff, £ = 4, . . . , 20 
to the Lense-Thirring combined precessions. 



i 


ACm (GGM02S-ITG-Gracc02s) 


ox + cry 


fi (mas yr" 1 ) 


4 


1.9 x 10 -11 


8.7 x 10~ 12 


7.2 


6 


2.1 x 10 -11 


4.6 x 10~ 12 


4.6 


8 


5.7 x 10~ 12 


2.8 x 10~ 12 


0.2 


10 


4.5 x 10~ 12 


2.0 x 10~ 12 




12 


1.5 x 10 -12 


1.8 x 10~ 12 




14 


6.6 x 10" 12 


1.6 x 10~ 12 




16 


2.9 x 10~ 12 


1.6 x 10~ 12 




18 


1.4 x 10~ 12 


1.6 x 10~ 12 




20 


2.0 x 10~ 12 


1.6 x 10 -12 






Sfi = 25% (SAV) 


Sfi = 18% (RSS) 





Table 3 Bias due to the mismodelling in the even zonals of the models X=ITG-Grace03s (Mayer-Gurr 2007), based on 
GRACE-only accumulated normal equations from data out of September 2002- April 2007 (neither apriori information nor 
regularization used), and Y=GGM02S (Tapley et al. 2005). The a for both models are formal. ACm are always larger 
than the linearly added sigmas, apart from I = 12 and £ — 18. 



i 


AC OT (ITG-Grace03s-GGM02S) 


ox + cry 


fi (mas yr" 1 ) 


4 


2.58 x 10 -11 


8.6 x 10~ 12 


9.6 


6 


1.39 x 10 -11 


4.7 x 10~ 12 


3.1 


8 


5.6 x 10~ 12 


2.9 x 10~ 12 


0.2 


10 


1.03 x 10 -11 


2 x 10" 12 




12 


7 x 10" 13 


1.8 x 10~ 12 




14 


7.3 x 10 -12 


1.6 x 10~ 12 




16 


2.6 x 10~ 12 


1.6 x 10~ 12 




18 


8 x 10~ 13 


1.6 x 10~ 12 




20 


2.4 x 10 -12 


1.6 x 10~ 12 






Sfi = 27% (SAV) 


<5a* = 21% (RSS) 





1.5 x 1CT 5 , as claimed by Ciufolini et al. (2009); thus, 
according to their reasoning corrected for this error, our 
approach would be exactly equivalent to state that the 
present-day accuracy in G would be wrong by a fac- 
tor of 100. Moreover, the relative uncertainties in the 
Earth's gravity models considered arc, in fact, all of the 
same order of magnitude; for example, the relative un- 
certainty in C 40 is 7.2 x 10~ 6 from EIGEN-GRACE02S 
(Rcigber et al. 2005) and 7.8 x 10 -6 in the more recent 
GGM03S model 21 (Tapley et al. 2007). Thus, the com- 
parison drawn by Ciufolini et al. (2009) between G and 
the even zonals is misleading. Even more important, it 
is well known that the rejection of a "suspect" exper- 



21 The relative uncertainty in C40 is 2 X 1CT 7 for ITG-Grace2010s 
(Mayer-Gurr et al. 2010), but it must be recalled that for such a 
model the available errors are the formal, statistical ones. 



imental result from a sample of data is always a very 
delicate matter (Taylor 1997), and quantitative criteria 
are needed to reject one or more outliers (Peirce 1852; 
Chauvenet 1863). Instead, Ciufolini et al. (2009) do not 
apply any of them to support their claims against the 
comparison of various geopotcntial models. The most 
famous rejection criterion is perhaps the one devised 
by Chauvenet (1863). Anyway, it relics upon an arbi- 
trary assumption that a measurement may be rejected 
if the probability of obtaining the deviation from the 
mean for that value is less than the inverse of twice the 
number of measurements; moreover, it makes no dis- 
tinction between the case of one or several suspicious 
data values. The criterion by Peirce (1852), instead, is 
a rigorous theory that can be easily applied in the case 
of several suspicious data values using the table in Ross 
(2003). Let us apply it to the case of G. If we consider 
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Table 4 Bias due to the mismodelling in the even zonals of the models X = GGM02S (Tapley et al. 2005) and Y = 
GGM03S (Tapley et al. 2007) retrieved from data spanning January 2003 to December 2006. The a for GGM03S are 
calibrated. ACeo are larger than the linearly added sigmas for £ = 4, 6. (The other zonals are of no concern) 



£ 


AC m (GGM02S-GGM03S) 


ox + ov 


ft (mas yr^ 1 ) 


4 


1.87 x 10 -11 


1.25 x 10 -11 


6.9 


6 


1.96 x 10 -11 


6.7 x 10~ 12 


4.2 


8 


3.8 x 10~ 12 


4.3 x 10~ 12 


0.1 


10 


8.9 x 10" 12 


2.8 x 10~ 12 


0.1 


12 


6 x 10" 13 


2.4 x 10 -12 




14 


6.6 x 10 -12 


2.1 x 10~ 12 




16 


2.1 x 10~ 12 


2.0 x 10~ 12 




18 


1.8 x 10~ 12 


2.0 x 10~ 12 




20 


2.2 x 10 -12 


1.9 x 10~ 12 






Sfi = 24% (SAV) 


Sfi = 17% (RSS) 





Table 5 Bias due to the mismodelling in the even zonals of the models X = EIGEN-GRACE02S (Reigber et al. 2005) 
and Y = GGM03S (Tapley et al. 2007). The o for both models are calibrated. AC to are always larger than the linearly 
added sigmas apart from I = 14, 18. 



t 


AC«) (EIGEN-GRACE02S-GGM03S) 


Ox + OY 


// (mas yr" 1 ) 


4 


2.00 x 10 -11 


8.1 x 10~ 12 


7.4 


6 


2.92 x 10 -11 


4.3 x 10~ 12 


6.3 


8 


1.05 x 10" 11 


3.0 x 10~ 12 


0.4 


10 


7.8 x 10~ 12 


2.9 x 10~ 12 


0.1 


12 


3.9 x 10~ 12 


1.8 x 10~ 12 




14 


5 x 10 -13 


1.7 x 10" 12 




16 


1.7 x 10~ 12 


1.4 x 10~ 12 




18 


2 x 10 -13 


1.4 x 10~ 12 




20 


2.5 x 10" 12 


1.4 x 10~ 12 






Sfi = 30% (SAV) 


Sfi = 20% (RSS) 





the set of modern measurements in Table 13 of Mohr 
and Taylor (1999), accurate to 10~ 3 - 1CT 4 , and the 
result by Cavendish (1798), accurate to 10 -2 , reported 
by dc Boer (1984) and Ohanian and Ruffini (1994), it 
turns out that, according to the Peirce (1852) criterion, 
the oldest value must be rejected. On the contrary, if 
we apply the Peirce (1852) criterion to all the values of, 
e.g., C40 from the models considered, no one of them 
has to be rejected. Interestingly, if we also considered 
models including CHAMP, LAGEOS and Earth-based 
data, and even the latest SLR-based solutions, the re- 
sult would not change: indeed, concerning C40, only the 
CHAMP-based TUM-2S model (Wermuth et al. 2004) 
would not pass the Peirce (1852) criterion. Thus, we 
do not see any founded, quantitative reasons to decline 
the comparison among different Earth's gravity models 
followed here. 

Another way to evaluate the uncertainty in the 
LAGEOS-LAGEOS II node test may consist of com- 
puting the nominal values of the total combined pre- 
cessions for different models and comparing them, i.e. 



by taking 

£ ^LAGEOS + Ci QLAGEOS II) ^ (X) _ Jg{Y)] 

(37) 

The results for each pair of models are shown in Table 
14. Their average is about 17%. 

A further, different approach that could be followed 
to take into account the scatter among the various solu- 
tions consists in computing mean and standard devia- 
tion of the entire set of values of the even zonals for the 
models considered so far, degree by degree, and taking 
the standard deviations as representative of the uncer- 
tainties S Jt, i = 4, 6, 8, ... . It yields Sfi = 15%, a figure 
slightly larger that that by Rics et al. (2008). Anyway, 
in evaluating mean and standard deviation for each 
even zonals, Ries et al. (2008) also used global grav- 
ity solutions like EIGEN-GL04C (Forste et al. 2006) 
and EIGEN-GL05C (Forste et al. 2008) which include 
data from the LAGEOS satellite itself; this may likely 
have introduced a sort of favorable a priori "imprint" 
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Table 6 Bias due to the mismodelling in the even zonals of the models X = JEM01-RL03B, based on 49 months of 
GRACE-only data, and Y = GGM03S (Tapley et al. 2007). The a for GGM03S are calibrated. AC m are always larger 
than the linearly added sigmas apart from £ = 16. 



1 


AC m (JEM01-RL03B-GGM03S) 


fx + cry 


ft (mas yr~ : ) 


4 


1.97 x 10 -11 


4.3 x 10~ 12 


7.3 


6 


2.7 x 10~ 12 


2.3 x 10~ 12 


0.6 


8 


1.7 x 10~ 12 


1.6 x 10~ 12 




10 


2.3 x 10~~ 12 


8 x 10" 13 




12 


7 x 10" 13 


7 x 10" 13 




14 


1.0 x 10 -12 


6 x 10~ 13 




16 


2 x 10~ 13 


5 x 10~ 13 




18 


7 x 10~ 13 


5 x 10~ 13 




20 


5 x 10" 13 


4 x 10" 13 






Sfi = 17% (SAV) 


5ft = 15% (RSS) 





Table 7 Bias due to the mismodelling in the even zonals of the models X = JEM01-RL03B and Y = ITG-Grace03s 
(Mayer-Giirr 2007). The a for ITG-Grace03s are formal. ACm are always larger than the linearly added sigmas. 



t 


AC OT (JEM01-RL03B-ITG-Grace03s) 


0"X + CTY 


ft (mas yr" 1 ) 


4 


2.68 x 10 -11 


4 x 10~ 13 


9.9 


6 


3.0 x 10~ 12 


2 x 10~ 13 


0.6 


8 


3.4 x 10~ 12 


1 x 10~ 13 


0.1 


10 


3.6 x 10 -12 


1 x 10" 13 




12 


6 x 10" 13 


9 x 10" 14 




14 


1.7 x 10~ 12 


9 x 10" 14 




16 


4 x 10~ 13 


8 x 10~ 14 




18 


4 x 10~ 13 


8 x 10~ 14 




20 


7 x 10" 13 


8 x 10" 14 






Sfi = 22% (SAV) 


Sfi = 10% (RSS) 





of the Lense-Thirring effect itself. Moreover, Ries et al. 
(2008) gave only a RSS evaluation of the total bias. 

It should be recalled that also the further bias due to 
the cross-coupling between J2 and the orbit inclination, 
evaluated to be about 9% in Iorio (2007d), must be 
added. 

2.1.2 An a-priori, "imprinting" effect? 

GRACE recovers the spherical harmonic coefficients of 
the gcopotcntial from the tracking of both satellites by 
GPS and from the observed intersatellite distance vari- 
ations (Reigber et al. 2005). A potential critical issue 
is, thus, a possible "memory" effect of the gravitomag- 
netic force. Its impact in the satellite-to-satellite track- 
ing was preliminarily addressed in Iorio (2005a); here 
we will focus on the "imprint" coming from the GRACE 
orbits which is more important for us because it mainly 
resides in the low degree even zonals (Iorio 2010b). 

Concerning that issue, Ciufolini and Pavlis (2005) 
write that such a kind of leakage of the Lense-Thirring 
signal itself into the even zonals retrieved by GRACE 
is completely negligible because the GRACE satellites 
move along (almost) polar orbits. Indeed, for perfectly 
polar (/ = 90 deg) trajectories, the gravitomagnctic 



force is entirely out-of-plane, while the perturbing ac- 
tion of the even zonals is confined to the orbital plane 
itself. According to Ciufolini and Pavlis (2005), the de- 
viations of the orbit of GRACE from the ideal polar or- 
bital configuration would have negligible consequences 
on the "imprint" issue. In particular, they write: "the 
values of the even zonal harmonics determined by the 
GRACE orbital perturbations are substantially inde- 
pendent on the a priori value of the LenseThirring ef- 
fect. [...] The small deviation from a polar orbit of 
the GRACE satellite, that is 1.7 x 10 -2 rad, gives only 
rise, at most, to a very small correlation with a fac- 
tor 1.7 x 10 -2 ". The meaning of such a statement is 
unclear; anyway, we will show below that such a con- 
clusion is incorrect. 

The relevant orbital parameters of GRACE are 
quoted in Table 1; variations of the orders of about 
10 km in the semimajor axis a and 0.001 deg in the 
inclination / may occur, but it turns out that they 
are irrelevant in our discussion (http://www.csr. utexas. 
edu/grace/ground/globe.html). The orbital plane of 
GRACE is, in fact, shifted by 0.98 deg from the ideal 
polar configuration, and, contrary to what is claimed in 
Ciufolini and Pavlis (2005), this docs matter because 
its classical secular node precessions are far from be- 
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Table 8 Aliasing effect of the mismodelling in the even zonal harmonics estimated in the X=ITG-Grace03s (Mayer-Giirr 
2007) and the Y=EIGEN-GRACE02S (Reigber et al. 2005) models. The covariance matrix a for ITG-Grace03s are formal, 
while the ones of EIGEN-GRACE02S are calibrated. AC to are larger than the linearly added sigmas for I = 4, 20, apart 
from I = 18. 



t 


ACffl (ITG-Grace03s-EIGEN-GRACE02S) 


fx + cr Y 


fe (mas yr" 1 ) 


4 


2.72 x 10" 11 


3.9 x 10" 12 


10.1 


6 


2.35 x 10" 11 


2.0 x 10~ 12 


5.1 


8 


1.23 x 10" 11 


1.5 x 10~ 12 


0.4 


10 


9.2 x 10~ 12 


2.1 x 10~ 12 


0.1 


12 


4.1 x 10" 12 


1.2 x 10~ 12 




14 


5.8 x 10" 12 


1.2 x 10~ 12 




16 


3.4 x 10" 12 


9 x 10 -13 




18 


5 x 10" 13 


1.0 x 10~ 12 




20 


1.8 x 10" 12 


1.1 x 10~ 12 






6{i = 37% (SAV) 


8[i = 24% (RSS) 





Table 9 Bias due to the mismodelling in the even zonals of the models X = JEM01-RL03B, based on 49 months of 
GRACE-only data, and Y = AIUB-GRACE01S (Jaggi et al. 2008). The latter one was obtained from GPS satellite-to- 
satellite tracking data and K-band range-rate data out of the period January 2003 to December 2003 using the Celestial 
Mechanics Approach. No accelerometer data, no de-aliasing products, and no regularisation was applied. The a for 
AIUB-GRACE01S are formal. AC to are always larger than the linearly added sigmas. 



1 


ACio (JEM01-RL03B— AIUB-GRACE01S) 


C"X + cry 


// (mas yr -1 ) 


4 


2.95 x 10" 11 


2.1 x 10~ 12 


11 


6 


3.5 x 10~ 12 


1.3 x 10~ 12 


0.8 


8 


2.14 x 10" 11 


5 x 10" 13 


0.7 


10 


4.8 x 10" 12 


5 x 10" 13 




12 


4.2 x 10~ 12 


5 x 10~ 13 




14 


3.6 x 10~ 12 


5 x 10~ 13 




16 


8 x 10~ 13 


5 x 10~ 13 




18 


7 x 10" 13 


5 x 10" 13 




20 


1.0 x 10" 12 


5 X 10" 13 






Sfj, = 26% (SAV) 


Sfi = 23% (RSS) 





ing negligible with respect to our issue. The impact of 
the Earth's gravitomagnetic force on the even zonals 
retrieved by GRACE can be quantitatively evaluated 
by computing the "effective" value C eo of the normal- 
ized even zonal gravity coefficients which would induce 
classical secular node precessions for GRACE as large 
as those due to its Lense-Thirring effect, which is in- 
dependent of the inclination /. To be more precise, 
C' m comes from solving the following equation which 
connects the classical even zonal precession of degree I 
ttj e = fl.eJi to the Lense-Thirring node precession IIlt 

ti.iJi = n LT . (38) 

Table 15 lists C eo for degrees I = 4, 6, which are the 
most effective in affecting the combination of eq. (32). 
Thus, the gravitomagnetic field of the Earth contributes 
to the value of the second even zonal of the geopoten- 
tial retrieved from the orbital motions of GRACE by 
an amount of the order of 2 x 10 -10 , while for I = 6 
the imprint is one order of magnitude smaller. Given 



the present-day level of accuracy of the latest GRACE- 
bascd solutions, which is of the order of 10 -12 , effects 
as large as those of Table 15 cannot be neglected. Thus, 
we conclude that the influence of the Earth's gravito- 
magnetic field on the low-degree even zonal harmonics 
of the global gravity solutions from GRACE may exist, 
falling well within the present-day level of measurabil- 
ity 

A further, crucial step consists of evaluating the im- 
pact of such an a-priori "imprint" on the test conducted 
with the LAGEOS satellites and the combination of 
eq. (32): if the LAGEOS-LAGEOS II uncancelled com- 
bined classical geopotential precession computed with 
the GRACE-based a-priori "imprinted" even zonals of 
Table 15 is a relevant part of - or even larger than - the 
combined Lense-Thirring precession, it will be demon- 
strated that the doubts concerning the a-priori gravit- 
omagnetic "memory" effect are founded. It turns out 
that this is just the case because eq. (32) and Table 15 
yield a combined geopotential precession whose magni- 
tude is 77.8 mas yr -1 (—82.9 mas yr -1 for I = 4 and 
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Table 10 Bias due to the mismodelling in the even zonals of the models X = EIGEN-GRACE02S (Reigber et al. 2005) 
and Y = AIUB-GRACE01S (J&ggi et al. 2008). The a for AIUB-GRACE01S are formal, while those of EIGEN-GRACE02S 
are calibrated. AC<o are larger than the linearly added sigmas for I = 4, 6, 8, 16. 



£ 


AC«) (EIGEN-GRACE02S— AIUB-GRACE01S) 


cx + OY 


ft (mas yr J ) 


4 


2.98 x 10 -11 


6.0 x 10~ 12 


11.1 


6 


2.29 x lO" 11 


3.3 x 10~ 12 


5.0 


8 


1.26 x 10" 11 


1.9 x 10~ 12 


0.4 


10 


6 x 10" 13 


2.5 x 10~ 12 




12 


5 x 10" 13 


1.6 x 10~ 12 




14 


5 x 10~ 13 


1.6 x 10~ 12 




16 


2.9 x 10~ 12 


1.4 x 10~ 12 




18 


6 x 10~ 13 


1.4 x 10~ 12 




20 


2 x 10" 13 


1.5 x 10" 12 






Sfi = 34% (SAV) 


Sfi = 25% (RSS) 





Table 11 Bias due to the mismodelling in the even zonals of the models X = JEM01-RL03B and Y = AIUB-GRACE02S 
(Jfiggi et al. 2009). The a for both AIUB-GRACE02S and JEM01-RL03B are formal. ACeo are larger than the linearly 
added sigmas for £ = 4 — 20. 



1 


AC eo (JEM01-RL03B— AIUB-GRACE02S) 


<Tx + cr Y 


ft (mas yr" ) 


4 


1.58 x 10 -11 


2 x 10" 13 


5.9 


6 


5.9 x 10" 12 


1 x 10~ 13 


1.3 


8 


5.8 x 10" 12 


1 x 10" 13 


0.2 


10 


1.27 x 10" 11 


1 x 10" 13 


0.1 


12 


4.3 x 10" 12 


1 x 10" 13 




14 


2.7 x 10~ 12 


1 x 10~ 13 




16 


1.6 x 10~ 12 


1 x 10~ 13 




18 


1.8 x 10~ 12 


1 x 10~ 13 




20 


1.9 x 10" 12 


1 x 10" 13 






Sfi = 16% (SAV) 


Sfi = 13% (RSS) 





5.1 mas yr -1 for I = 6), i.e. just 1.6 times the Lense- 
Thirring signal itself. This means that the part of the 
LAGEOS-LAGEOS II uncancelled classical combined 
node precessions which is affected by the "imprinting" 
of the Lense-Thirring force through the GRACE-based 
geopotential's spherical harmonics is as large as the 
LAGEOS-LAGEOS II combined gravitomagnetic sig- 
nal itself. 

We, now, comment on how Ciufolini and Pavlis 
(2005) reach a different conclusion. They write: "How- 
ever, the Lense-Thirring effect depends on the third 
power of the inverse of the distance from the central 
body, i.e., (1/r) 3 , and the J2, J4, Jq,. . . effects depend 
on the powers (1/r) 3 - 5 , (1/r) 5 5 , (1/r) 7 5 ,. . . of the 
distance; then, since the ratio of the semimajor axes 
of the GRACE satellites to the LAGEOS' satellites is 
~ 1227V — 1-8, any conceivable "Lense-Thirring im- 
print" on the spherical harmonics at the GRACE alti- 
tude becomes quickly, with increasing distance, a neg- 
ligible effect, especially for higher harmonics of degree 
I > 4. Therefore, any conceivable "Lense-Thirring im- 
print" is negligible at the LAGEOS' satellites altitude." 
From such statements it seems that they compare the 
classical GRACE precessions with the gravitomagnetic 



LAGEOS' ones. This is meaningless since, as we have 
shown, one has, first, to compare the classical and rel- 
ativistic precessions of GRACE itself, with which the 
Earth's gravity field is solved for and only after to com- 
pute the impact of the relativistically "imprinted" part 
of the GRACE-based even zonals on the combined LA- 
GEOS nodes. These two stages have to be kept sepa- 
rate, with the first one which is fundamental; if different 
satellite(s) Y were to be used to measure the gravito- 
magnetic field of the Earth, the impact of the Lense- 
Thirring effect itself on them should be evaluated by 
using the "imprinted" even zonals evaluated in the first 
stage. Finally, in their latest statement Ciufolini and 
Pavlis (2005) write: "In addition, in (Ciufolini et al. 
1997b), it was proved with several simulations that by 
far the largest part of this "imprint" effect is absorbed 
in the by far largest coefficient J2." Also such a state- 
ment, in the present context, has no validity since the 
cited work refers to a pre-GRACE era. Moreover, no 
quantitative details at all were explicitly released con- 
cerning the quoted simulations, so that it is impossible 
to form an opinion. 
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Table 12 Bias due to the mismodclling in the even zonals of the models X = ITG-Grace2010s (Mayer-Giirr et al. 2010) and 
Y = AIUB-GRACE02S (Jaggi et al. 2009). The ITG-Grace2010s model has been obtained by processing 7 yr (2002-2009) 
of GRACE data. The a for both models are formal. ACeo are larger than the linearly added sigmas for all the degrees 



considered. 


1 AC m (ITG-Grace2010s-AIUB-GRACE02S) 


o"x + cry 


ft (mas yr - ) 


4 2.665 x 10 _li 


1.9 x 10~ la 


9.9 


6 4.66 x 10~ 12 


1.1 x 10~ 13 


1.0 


8 3.86 x 10~ 12 


9 x 10~ 14 


0.1 


10 9.50 x 10~ 12 


8 x 10~ 14 


0.1 


12 1.67 x 10" 12 


7 x RT 14 




14 2.36 x 10" 12 


7 x RT 14 




16 8.7 x 10~ 13 


6 x 10~ 14 




18 8.0 x 10" 13 


6 x 10~ 14 




20 1.09 x 10~ 12 


7 x RT 14 




Sfi = 23% (SAV) 


5fi = 21% (RSS) 





Table 13 Bias due to the mismodelling in the even zonals of the models X = ITG-Grace2010s (Mayer-Giirr et al. 2010) 
and Y = GGM03S (Tapley et al. 2007). The a of ITG-Grace2010s are formal, while those for GGM03S are calibrated. 
ACeo are larger than the linearly added sigmas for £ = 4, 10, 12, 16, 18. 



£ 


AC l0 (ITG-Grace2010s-GGM03S) 


o"x + cry 


ft (mas yr _I ) 


4 


3.05 x lO" 11 


4.3 x 10~ 12 


11.3 


6 


1.4 x 10~ 12 


2.3 x 10~ 12 


0.3 


8 


3 x 10~ 13 


1.6 x 10~ 12 




10 


9 x 10" 13 


8 x 10~ 13 




12 


1.9 x 10~ 12 


6 x 10~ 13 




14 


6 x 10" 13 


6 x 10~ 13 




16 


9 x 10 -13 


4 x 10~ 13 




18 


1.7 x 10 -12 


4 x 10~ 13 




20 


2 x 10" 13 


4 x 10" 13 






Sfi = 25% (SAV) 


Sfi = 24% (RSS) 





2.1.3 A new approach to extract the Lense-Thirring 
signature from the data 

The technique adopted so far by Ciufolini and Pavlis 
(2004) and Ries et al. (2008) to extract the gravito- 
magnetic signal from the LAGEOS and LAGEOS II 
data is described in detail in (Lucchesi and Balmino 
2006; Lucchesi 2007b). In both approaches the Lense- 
Thirring force is not included in the dynamical force 
models used to fit the satellites' data. In the data re- 
duction process no dedicated gravitomagnctic parame- 
ter is estimated, contrary to, e.g., station coordinates, 
state vector, satellite drag and radiation coefficients Co 
and Cr, respectively, etc.; its effect is retrieved with a 
sort of post-post-fit analysis in which the time series of 
the computed 22 "residuals" of the nodes with the differ- 
ence between the orbital elements of consecutive arcs, 
combined with eq. (32), is fitted with a straight line. 



The expression "residuals of the nodes" is used, strictly speak- 
ing, in an improper sense because the Keplerian orbital elements 
are not directly measured. 



In regard to possible other approaches which could 
be followed, it would be useful to, e.g., estimate (in 
the least square sense), among other solve- for param- 
eters, purely phenomenological corrections Af2 to the 
LAGEOS/LAGEOS II node precessions as well, with- 
out modelling the Lense-Thirring effect itself, so that 
it will be, in principle, contained in AO, and combine 
them according to eq. (32). Something similar has 
been done - although for different scopes - for the per- 
ihelia of the inner planets of the solar system (Pitjcva 
2005a) (see Section 3) and the periastron of the pul- 
sars (Kramer et al. 2006). To be more definite, various 
solutions with a complete suite of dynamical models, 
apart from the gravitomagnctic force itself, should be 
produced in which one inserts a further solve-for param- 
eter, i.e. a correction Af2 to the standard Newtonian 
modelled precessions. One could see how the outcome 
varies by changing the data sets and/or the parameters 
to be solved for. Maybe it could be done for each arc, 
so to have a collection of such node extra-rates. Such a 
strategy would be much more model-independent. 

As previously suggested by Nordtvedt (2001), an- 
other way to tackle the problem consists of looking at 
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Table 14 Systematic uncertainty S/j, in the LAGEOS-LAGEOS II test evaluated by taking the absolute value of the 
difference between the nominal values of the total combined node precessions due to the even zonals for different models X 
and Y, i.e. I fF copot (X) - Si gcopot (Y) I . The average is « 17%. 



Models compared S/i 

ITG-Grace2010s-JEM01-RL03B 7% 

ITG-Grace2010s-GGM02S 0.3% 

ITG-Grace2010s-GGM03S 24% 

ITG-Grace2010s-ITG-Grace02 25% 

ITG-Grace2010s-ITG-Grace03 27% 

ITG-Grace2010s-AIUB-GRACE02S 23% 

ITG-Grace2010s-AIUB-GRACE01S 27% 

ITG-Grace2010s-EIGEN-GRACE02S 5% 

AIUB-GRACE02S— JEM01-RL03B 16% 

AIUB-GRACE02S— GGM02S 23% 

AIUB-GRACE02S— GGM03S 17% 

AIUB-GRACE02S-ITG-Grace02 2% 

AIUB-GRACE02S-ITG-Grace03 13% 
AIUB-GRACE02S— EIGEN-GRACE02S 28% 

AIUB-GRACE01S— JEM01-RL03B 20% 

AIUB-GRACE01S— GGM02S 27% 

AIUB-GRACE01S— GGM03S 3% 

AIUB-GRACE01S-ITG-Grace02 2% 

AIUB-GRACE01S-ITG-Grace03 0.1% 
AIUB-GRACE01S— EIGEN-GRACE02S 33% 

JEM01-RL03B— GGM02S 7% 

JEM01-RL03B— GGM03S 17% 

JEM01-RL03B-ITG-Grace02 18% 

JEM01-RL03B-ITG-Grace03s 20% 

JEM01-RL03B— EIGEN-GRACE02S 13% 

GGM02S— GGM03S 24% 

GGM02S-ITG-Grace02 25% 

GGM02S-ITG-Grace03s 27% 

GGM02S— EIGEN-GRACE02S 6% 

GGM03S-ITG-Grace02 1% 

GGM03S-ITG-Grace03s 3% 

GGM03S— EIGEN-GRACE02S 30% 

ITG-Grace02-ITG-Grace03s 2% 

ITG-Gracc02-EIGEN-GRACE02S 31% 

ITG-Grace03s— EIGEN-GRACE02S 33% 



a Lense-Thirring-dcdicated parameter to be estimated 
along with all the zonals in a new global solution for the 
gravity field incorporating the gravitomagnctic compo- 
nent as well; instead, in all the so far produced global 
gravity solutions no relativistic parameter (s) have been 
included in the set of the estimated ones. As shown in 
Section 2.1.2, this would also cure the impact of possi- 
ble forms of a-priori "imprinting" effects. 

A first, tentative step towards the implementation 
of the strategy of the first point mentioned above with 
the LAGEOS satellites in terms of the PPN parameter 
7 has been recently taken by Combrinck (2008). 



2.2 The LARES mission 

Van Patten and Everitt (1976a) proposed to measure 
the Lcnsc-Thirring precession of the nodes f2 of a pair of 
counter-orbiting spacecraft to be launched in terrestrial 
polar orbits and endowed with a drag-free apparatus. 
A somewhat equivalent, cheaper version of such an idea 
was put forth ten years later by Ciufolini (1986, 1989) 
who proposed to launch a passive, geodetic satellite in 
an orbit identical to that of LAGEOS apart from the 
orbital planes which should have been displaced by 180 
deg apart. The measurable quantity was, in the case of 
the proposal by Ciufolini (1986), the sum of the nodes 
of LAGEOS and of the new spacecraft, later named 
LAGEOS III (Ciufolini 1994), LARES (Ciufolini et al. 
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Table 15 c\ : effective "gravitomagnetic" normalized gravity coefficients for GRACE (£ = 4, 6; m = 0). They have been 
obtained by comparing the GRACE classical node precessions to the Lense-Thirring rate. Thus, they may be viewed as a 
quantitative measure of the leakage of the Lense-Thirring effect itself into the second and third even zonal harmonics of the 
global gravity solutions from GRACE. Compare them with the much smaller calibrated errors in C40 and C§o of the 
GGM03S model (Tapley et al. 2007). 

2.23 x irr 10 -2.3 x lO" 11 4 x 10~ 12 2 x lfr 12 



1998b), and WEBER-SAT 23 , in order to cancel to a 
high level of accuracy the corrupting effect of the mul- 
tipoles of the Newtonian part of the terrestrial gravi- 
tational potential which represent the major source of 
systematic errors. Although extensively studied by var- 
ious groups (Ries et al. 1989; Ciufolini et al. 1998b), 
such an idea was not implemented for many years. Io- 
rio et al. (2002a) proposed to include also the data from 
LAGEOS II by using a different observable. Such an 
approach was proven in Iorio (2005b) to be, in prin- 
ciple, potentially useful in making the constraints on 
the orbital configuration of the new SLR satellite less 
stringent than it was originally required in view of the 
recent improvements in our knowledge of the classical 
part of the terrestrial gravitational potential due to the 
dedicated CHAMP and, especially, GRACE missions. 

Since reaching high altitudes and minimizing the un- 
avoidable orbital injection errors is expensive, the pos- 
sibility of discarding LAGEOS and LAGEOS II using 
a low-altitude, nearly polar orbit for LARES (Lucchesi 
and Paolozzi 2001; Ciufolini 2006) was explored. How- 
ever, in Iorio (2002, 2007b) it was proven that such 
alternative approaches are not feasible. It was also sug- 
gested that LARES would be able to probe alternative 
theories of gravity (Ciufolini 2004), but also in this case 
it turned out to be impossible (Iorio 2005c, 2007c). 

The stalemate came to an end when ASI recently 
approved the LARES mission, although with a dif- 
ferent orbital geometry with respect to the original 
configuration: now the orbital altitude is 1450 km 
corresponding to a semimajor axis a = 7828 km 
(Ciufolini et al. 2009). See Table 1 for the new 
orbital parameters and the related Lense-Thirring 
node precession. LARES should be launched in late 
2010/carly 2011 with the first qualification flight of 
the VEGA rocket (http://www.spacenews.com/civil/ 
100115- asi- expects- budget- remain- flat- 20 1 . html) . 

The combination that should be used for measuring 
the Lense-Thirring effect with LAGEOS, LAGEOS II 

23 In memory of Dr J. Weber, US Naval Academy (USNA) Class 
of 1940. 



and LARES is (Iorio 2005b) 

^LAGEOS + ^ ^LAGEOS II + fc^LARES^ (39) 

where the coefficients k\ and ki entering eq. (39) are 
defined as 



fci 



LARES j^LAGEOS f|LAGEOS qLARES 



JJLAGEOS IIqLARES _^LARES ^LAGEOS II 



0.3586, 



^LAGEOS qLAGEOS II_^LAGEOS II^LAGEOS 
qLAGEOS II qLAHES _ nLARES nLACiEOS II 



0.0751. 



(40) 



By construction, the combination eq. (39) cancels the 
impact of the first two even zonals; we have used o,lr = 
7828 km and 7 LR = 71.5 deg. The total Lense-Thirring 
effect, according to eq. (39) and eq. (40), amounts to 
50.8 mas yr _1 . 

2.2.1 A conservative evaluation of the impact of the 
geopotential on the LARES mission 

The systematic error due to the uncancelled even zonals 
Jq^Jst-- can be conservatively evaluated as 



l^ 1 



LAGEOS + ftLAGEOS II + ^qLARES 



8J t 
(41) 



Of crucial importance is how to assess SJi. By pro- 
ceeding as in Section 2.1.1 and by using the same mod- 
els up to degree I = 60 because of the lower altitude 
of LARES with respect to LAGEOS and LAGEOS II 
which brings into play more even zonals, we come up 
with the results presented in Table 16. They have been 
obtained with the standard and widely used Kaula ap- 
proach (Kaula 1966) in the following way. We, first, 
calibrated our numerical calculation with the analyt- 
ical ones performed with the explicit expressions for 
Clj worked out up to £ = 20 in Iorio (2003c); then, 
after having obtained identical results, we confidently 
extended our numerical calculation to higher degrees by 
means of two different softwares (Matlab and MATH- 
EMATICAL 
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Table 16 Systematic percent uncertainty Sfi in the combined Lense-Thirring effect with LAGEOS, LAGEOS II and 
LARES according to eq. (41) and SJi = AJe up to degree £ — 60 for the global Earth's gravity solutions considered here; 
the approach by (Kaula 1966) has been followed. For LARES we adopted olr = 7828 km, Tlr = 71.5 deg, eLR = 0.0. 



Models compared (SJ e = AJ e ) 5ji (SAV) Sfj, (RSS) 



ITG-Grace2010s-JEM01-RL03B 


4% 


2% 


ITG-Grace2010s-GGM02S 


14% 


8% 


ITG-Grace2010s-GGM03S 


2% 


1% 


ITG-Grace2010s-ITG-Grace02 


■2% 


1% 


ITG-Grace2010s-ITG-Grace03 


2% 


1% 


ITG-Grace2010s-AIUB-GRACE02S 


5% 


2% 


ITG-Grace2010s-AIUB-GRACE01S 


22% 


14% 


ITG-Grace2010s-EIGEN-GRACE02S 


24% 


13% 


AIUB-GRACE02S— JEM01-RL03B 


9% 


4% 


AIUB-GRACE02S— GGM02S 


16% 


7% 


AIUB-GRACE02S— GGM03S 


4% 


2% 


AIUB-GRACE02S-ITG-Grace02 


5% 


2% 


AIUB-GRACE02S-ITG-Grace03 


6% 


3% 


AIUB-GRACE02S— EIGEN-GRACE02S 


24% 


13% 


AIUB-GRACE02S— AIUB-GRACE01S 


23% 


13% 


AIUB-GRACE01S— JEM01-RL03B 


23% 


16% 


AIUB-GRACE01S— GGM02S 


16% 


8% 


AIUB-GRACE01S— GGM03S 


22% 


13% 


AIUB-GRACE01S-ITG-Grace02 


24% 


15% 


AIUB-GRACE01S-ITG-Grace03 


22% 


14% 


AIUB-GRACE01S— EIGEN-GRACE02S 


14% 


7% 


JEM01-RL03B— GGM02S 


14% 


9% 


JEM01-RL03B— GGM03S 


5% 


3% 


JEM01-RL03B-ITG-Grace02 


4% 


2% 


JEM01-RL03B-ITG-Grace03s 


5% 


2% 


JEM01-RL03B— EIGEN-GRACE02S 


26% 


15% 


GGM02S— GGM03S 


13% 


7% 


GGM02S-ITG-Grace02 


16% 


8% 


GGM02S-ITG-Grace03s 


14% 


7% 


GGM02S— EIGEN-GRACE02S 


14% 


7% 


GGM03S-ITG-Grace02 


3% 


2% 


GGM03S-ITG-Grace03s 


2% 


0.5% 


GGM03S— EIGEN-GRACE02S 


24% 


13% 


ITG-Grace02-ITG-Grace03s 


3% 


2% 


ITG-Grace02— EIGEN-GRACE02S 


25% 


14% 


ITG-Grace03s— EIGEN-GRACE02S 


24% 


13% 



It must be stressed that our results may be still op- 
timistic: indeed, computations for £ > 60 become un- 
reliable because of numerical instability of the results. 

In Table 17 we repeat the calculation by using for SJg 
the covariance matrix sigmas <jj t ; also in this case we 
use the approach by Kaula (1966) up to degree I = 60. 

If, instead, one assumes SJ( = at, £ = 2,4,6, ... 
i.e., the standard deviations of the sets of all the best 
estimates of Ji for the models considered here the sys- 
tematic bias, up to £ = 60, amounts to 12% (SAV) and 
6% (RSS). Again, also this result may turn out to be 
optimistic for the same reasons as before. 

It must be pointed out that the evaluations pre- 
sented here rely upon calculations of the coefficients 
f2 j performed with the well known standard approach 



by Kaula (1966); it would be important to try to fol- 
low also different computational strategies in order to 
test them. In regard to this point, Ciufolini et al. (2009) 
state that, in reality, the bias due to the even zonals is of 
the order of 1% or less, and repeatedly write that the re- 
sults of one of us (L. I.) arc based on "miscalculations". 
In fact, Ciufolini ct al. (2009), do not demonstrate their 
allegations by explicitly disclosing the alleged error(s). 

2.2.2 The impact of some non- gravitational 
perturbations 

It is worthwhile noting that also the impact of the 
subtle non-gravitational perturbations will be different 
with respect to the original proposal because LARES 
will fly in a lower orbit and its thermal behavior will 
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Table 17 Systematic percent uncertainty Sfi in the combined Lense-Thirring effect with LAGEOS, LAGEOS II and 
LARES according to eq. (41) and SJe = <jj e up to degree I = 60 for the global Earth's gravity solutions considered here; 
the approach by (Kaula 1966) has been followed. For LARES we adopted olr = 7828 km, Tlr = 71.5 deg, clr = 0.0. 



Model (SJ e = a e ) Sfi (SAV) Sjj, (RSS) 



ITG-Grace2010s (formal) 


0.2% 


0.1% 


AIUB-GRACE02S (formal) 


1% 


0.9% 


AIUB-GRACE01S (formal) 


11% 


9% 


JEM01-RL03B (formal) 


1% 


0.9% 


GGM03S (calibrated) 


5% 


4% 


GGM02S (formal) 


20% 


15% 


ITG-Grace03s (formal) 


0.3% 


0.2% 


ITG-Grace02s (formal) 


0.4% 


0.2% 


EIGEN-GRACE02S (calibrated) 


21% 


17% 



probably be different with respect to LAGEOS and LA- 
GEOS II. The reduction of the impact of the thermal 
accelerations, like the Yarkovsky-Schach effects, should 
have been reached with two concentric spheres. How- 
ever, as explained by Andres (2007), this solution will 
increase the floating potential of LARES because of the 
much higher electrical resistivity and, thus, the pertur- 
bative effects produced by the charged particle drag. 
Moreover, the atmospheric drag will increase also be- 
cause of the lower orbit of the satellite, both in its neu- 
tral and charged components. Indeed, although it does 
not affect directly the node il, it induces a secular de- 
crease of the inclination / of a LAGEOS-likc satellite 
(Milani ct al. 1987) which translates into a further bias 
for the node itself according to 

in which 51 accounts not only for the measurement 
errors in the inclination, but also for any unmod- 
elled/mismodclled dynamical effect on it. According 
to (Iorio 2010c), the secular decrease for LARES would 
amount to 

(^-) » -0.6 mas yr" 1 (43) 

\ dt /LR 

yielding a systematic uncertainty in the Lense-Thirring 
signal of eq. (39) of about 3 — 9% yr _1 . An analogous 
indirect node effect via the inclination could be induced 
by the thermal Yarkovski-Rubincam force as well (Iorio 
2010c). Also the Earth's albedo, with its anisotropic 
components, may have a non-negligible effect. 

Ciufolini et al. (2009) objected that, in fact, the dis- 
turbing effect examined would not appear in the real 
data analysis procedure because the inclination along 
with all the other Kcplerian orbital elements would be 
measured arc by arc, so that one should only have to 
correct the signal for the measured value of the incli- 
nation; after all, the same problems, if not even larger, 



would occur with the semimajor axes of the LAGEOS 
satellites, which are known to undergo still unexplained 
secular decrease of 1.1 mm d _1 (Rubincam 1982) and 
their consequent mappings onto the node rates. The 
problem is that while a perturbation Aa pertains the in- 
plane, radial component (Christodoulidis et al. 1988) of 
the LAGEOS orbits, both the Lense-Thirring node pre- 
cession and the shifts in the inclination affect the out-of- 
plane, normal component of the orbit (Christodoulidis 
et al. 1988); thus, even if repeated corrections to the 
semimajor axis could be applied without affecting the 
gravitomagnetic signal of interest, the same would not 
hold for the inclination. This is particularly true in view 
of the fact that, for still unexplained reasons, the Lense- 
Thirring effect itself has never been estimated, either as 
a short-arc or as a global parameter. Moreover, it has 
been claimed that the recent improvements in atmo- 
spheric refraction modelling would allow to measure the 
inclination of the LAGEOS satellites at a level of accu- 
racy, on average, of 30 fi&s for LAGEOS and 10 /ias for 
LAGEOS II (Ciufolini ct al. 2009). Firstly, the track- 
ing of a relatively low satellite is always more difficult 
than for higher targets, so that caution would be needed 
in straightforwardly extrapolating results valid for LA- 
GEOS to the still non-existing LARES. Second, it is 
difficult to understand the exact sense of such claims 
because they would imply an accuracy 6r ss aSI in re- 
constructing the orbits of LAGEOS and LAGEOS II, 
on average, of 0.2 cm and 0.06 cm, respectively. 

Let us point out the following issue as well. At 
present, it is not yet clear how the data of LAGEOS, 
LAGEOS II and LARES will be finally used by the 
proponent team in order to detect the Lense-Thirring 
effect. This could turn out to be a non-trivial matter 
because of the non-gravitational perturbations. Indeed, 
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if, for instance, a combination 24 



^LARES + ^LAGEOS + ^LAGEOS II 



(44) 



was adopted instead of that of eq. (39), the coefficients 
of the nodes of LAGEOS and LAGEOS II, in view of 
the lower altitude of LARES, would be 



hi 



n 



LAGEOS II^LARES oLARESoLAGEOS II 



<1 L . 



^LARES j'jLAGEOS II_JjLAGEOS II^LAGEOS 



LARES ^LAGEOS oLAGEOS qLARES 



^LAGEOS^LAGEOS II_QLAGEOS II^LAGEOS 



13.3215, 



4.7744. 



(45) 



and the combined Lense-Thirring signal would amount 
to 676.8 mas yr _1 . As a consequence, the direct and 
indirect effects of the non-gravitational 25 perturbations 
on the nodes of LAGEOS and LAGEOS II would be 
enhanced by such larger coefficients and this may yield 
a lower total obtainable accuracy. 



3 In search of the Sun's gravitomagnetic field 

3.1 General considerations 

Recent determinations of the Sun's proper angular mo- 
mentum 



70 



(190.0 ± 1.5) x lO^ 9 kg in s 



(46) 



from helioseismology (Pijpcrs 1998, 2003), accurate 
to 0.8%, yield a value about one order of magnitude 
smaller than that obtained by assuming a homogeneous 
and uniformly rotating Sun, as done in the pioneer- 
ing work by de Sitter (1916), and also by Soffel (1989) 
and Cugusi and Proverbio (1978) who concluded that, 
at their time, it was not possible to measure the solar 
Lense-Thirring effect. Despite the reduced magnitude 
of the solar gravitomagnetic field with respect to the 
earlier predictions, the present and near future situa- 
tion seems more promising. The characteristic length 
with which the accuracy of the determination of the 
orbits of the particles should be compared is 



l Q 



M Pi c 



= 319 m; 



(47) 



in the case of gravitoelectric effects, l g is usually re- 
placed by the Schwarzschild radius R g = 2GM/c 2 = 3 



km for the Sun. The present-day accuracy in knowing, 
e.g., the inner planets' mean orbital radius 



(r) = a [ 1 + - 



(48) 



is shown in Table 18. Such values have been obtained 
by linearly propagating the formal, statistical errors in 
a and e according to Table 3 of Pitjeva (2008); even by 
re-scaling them by a a factor of, say, 2 — 5, the grav- 
itomagnetic effects due to the Sun's rotation fall, in 
principle, within the measurability domain. Another 
possible way to evaluate the present-day uncertainty in 
the planetary orbital motions consists of looking at dif- 
ferent ephemcrides of comparable accuracy. In Table 18 
we do that for the EPM2006/EPM2008 (Pitjeva 2008, 
2010), and the DE414/DE421 (Standish 2006; Folkner 
et al. 2008) ephemcrides; although, larger than S(r), 
the maximum differences between such ephemcrides are 
smaller than the solar gravitomagnetic length I®. 

With regard to the currently collected ranging data 
to the Venus Express spacecraft 26 and in view of the 
ongoing Messenger (Balogh et al. 2007) and the future 
BepiColombo 27 (Milani et al. 2002; Balogh et al. 2007) 
missions to Mercury and of the developments in the 
Planetary Laser Ranging (PLR) technique, the plan- 
etary orbit accuracy is likely to be further improved. 
More specifically, recent years have seen increasing ef- 
forts towards the implementation of PLR accurate to 
cm-level (Smith et al. 2006; Chandler et al. 2004; Neu- 
mann et al. 2006; Degnan 2006; Turyshcv and Williams 
2007; Merkowitz et al. 2007; Degnan 2008; Zubcr and 
Smith 2008). It would allow to reach major improve- 
ments in three related fields: solar system dynamics, 
tests of general relativity and alternative theories of 
gravity, and physical properties of the target planet it- 
self. In principle, any solar system body endowed with 
a solid surface and a transparent atmosphere would be 
a suitable platform for a PLR system, but some targets 
are more accessible than others. Major efforts have 
been practically devoted so far to Mercury (Smith et 
al. 2006) and Mars (Chandler et al. 2004; Turyshev and 
Williams 2007), although simulations reaching 93 AU 
or more have been undertaken as well (Degnan 2006, 
2008). In 2005 two interplanetary laser transponder ex- 
periments were successfully demonstrated by the God- 
dard Geophysical Astronomical Observatory (GGAO). 
The first utilized the non-optimized Mercury Laser Al- 



24 The impact of the geopotential is, by construction, unaffected 
with respect to the combination of eq. (39) . 

25 The same may hold also for time-dependent gravitational per- 
turbations affecting the nodes of LAGEOS and LAGEOS II, like 
the tides. 



26 See on the WEB http://www.esa.int/esaMI/Venus_Express/. 

27 It is an ESA mission, including two spacecraft, one of which 
provided by Japan, to be put into orbit around Mercury. The 
launch is scheduled for 2014. The construction of the instruments 
is currently ongoing. 



22 



Table 18 First line: uncertainties (in m) in the average heliocentric distances of the inner planets obtained by propagating 
the formal errors in a and e according to Table 3 of Pitjeva (2008); the EPM2006 ephemerides were used by Pitjeva (2008). 
Second line: maximum differences (in m) between the EPM2006 and the DE414 (Standish 2006) ephemerides for the inner 
planets in the time interval 1960-2020 according to Table 5 of Pitjeva (2010). Third line: maximum differences (in m) 
between the EPM2008 (Pitjeva 2010) and the DE421 (Folkner et al. 2008) ephemerides for the inner planets in the time 
interval 1950-2050 according to Table 5 of Pitjeva (2010). They have to be compared with the characteristic gravitomagnetic 
length of the Sun if = 319 m. 



Type of orbit uncertainty 


Mercury 


Venus 


Earth 


Mars 


S(r) (EPM2006) 


38 


3 


1 


2 


EPM2006-DE414 


256 


131 


17.2 


78.7 


EPM2008-DE421 


185 


4.6 


11.9 
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timeter (MLA) on the Messenger spacecraft (Smith et 
al. 2006; Neumann et al. 2006), obtaining a formal er- 
ror in the laser range solution of 0.2 m, or one part 
in 10 11 . The second utilized the Mars Orbiting Laser 
Altimeter (MOLA) on the Mars Global Surveyor space- 
craft (Abshire 2006; Neumann et al. 2006). A precise 
value of the Earth-Mars distance, measured between 
their centers of mass and taken over an extended pe- 
riod (five years or more), would support, among other 
things, a better determination of several parameters of 
the solar system. Sensitivity analyses point towards 
measurement uncertainties between 1 mm and 100 mm 
(Chandler et al. 2004). A recent analysis of the fu- 
ture BepiColombo mission to Mercury, aimed to ac- 
curately determining, among other things, several key 
parameters of post-Newtonian gravity and the solar 
quadrupole moment from Earth-Mercury distance data 
collected with a multi-frequency radio link (Milani et 
al. 2002, 2008), points toward a maximum uncertainty 
of 4.5 — 10 cm in determining the Earth-Mercury range 
over a multi-year time span (1-8 yr) (Milani et al. 2002; 
Ashby et al. 2007; Milani et al. 2008). A proposed 
spacecraft-based mission aimed to accurately measure 
also the gravitomagnetic field of the Sun and its Ji 
along with other PPN parameters like 7 and j3 by means 
of interplanetary ranging is the Astrodynamical Space 
Test of Relativity using Optical Devices 28 (ASTROD) 
(Ni 2008). Another space- based mission proposed to 
accurately test several aspects of the gravitational in- 
teraction via interplanetary laser ranging is the Laser 
Astrometric Test of Relativity (LATOR) (Turyshev et 
al. 2009). 

It is remarkable to note that the currently available 
estimate of S@ from helioseismology is accurate enough 
to allow, in principle, a genuine Lense-Thirring test. 
Moreover, it was determined in a relativity-free fashion 
from astrophysical techniques which do not rely on the 



28 Its cheaper version ASTROD I makes use of one spacecraft in a 
Venus-gravity-assisted solar orbit, ranging optically with ground 
stations (Appourchaux et al. 2009). 



dynamics of planets in the gravitational field of the Sun. 
Thus, there is no any a priori "memory" effect of general 
relativity itself in the adopted value of Sq. 

3.2 The Lense-Thirring perihelion precessions of the 
planets and their mcasurability 

The action of the solar gravitomagnetic field on Mer- 
cury's longitude of perihelion 29 w was calculated for 
the first time by de Sitter (1916) who, by assuming a 
homogenous and uniformly rotating Sun, found a sec- 
ular rate of —0.01 arcseconds per century (arcsec cy" 1 
in the following). This value is also quoted by Soffcl 
(1989); Cugusi and Provcrbio (1978) yield -0.02 arcsec 
cy _1 for the argument of perihelion uj of Mercury. 

The recent estimate of the Sun's angular momentum 
(Pijpers 1998, 2003) from helioseismology implies a pro- 
cessional effect for Mercury which is one order of mag- 
nitude smaller; see Table 19 for the predicted Lense- 
Thirring precessions hjlt of the longitudes of the per- 
ihelia of the inner planets which are of the order of 
10~ 3 — 10~ 5 arcsec cy _1 . They are obtained by taking 
into account the fact that the inclinations / of the plan- 
ets usually quoted in the literature refer to the mean 
ecliptic at a given epoch 30 (Roy 2005), while the Sun's 
equator is tilted by e© = 7.15 deg to the mean ecliptic 
at J2000 (Beck and Giles 2005). So far, the solar Lense- 
Thirring effect on the orbits of the inner planets was 
believed to be too small to be detected (Soffcl 1989). 
However, the situation seems now favorably changing. 
Iorio (2005d) preliminarily investigated the possibility 
of measuring such tiny effects in view of recent impor- 
tant developments in the planetary ephemerides gener- 
ation. 

First attempts to measure the Sun's Lense-Thirring 
effect have recently been implemented by Iorio (2007e, 
2008a) with the Ephemerides of Planets and the Moon 



29 Since f2 and oj do not lie, in general, in the same plane, zu is a 
"dogleg" angle. 

30 It is J2000 (JD 2451545.0). 
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Table 19 First line: corrections Azu, in 10~ 4 arcsec cy -1 (1 arcsec cy -1 = 10 mas yr _1 = 1.5 x 1CP 15 s _1 ), to the 
standard Newton/Einstein secular precessions of the perihelia of the inner planets estimated by Pitjeva (2005a) with the 
EPM2004 ephemerides. The result for Venus (E.V. Pitjeva, private communication, 2008) has been obtained by recently 
processing radiometric data from Magellan spacecraft with the EPM2006 ephemerides (Pitjeva 2008); it has also been 
cited in Table 4 of Fienga et al. (2010) and, as far as the uncertainty is concerned, in Table 1 of Biswas and Mani (2008). 
The errors in square brackets are the 1 — a formal, statistical ones. Second line: predicted Lense-Thirring perihelion 
precessions tolt, in 10~ 4 arcsec cy~ x (Iorio 2007e). Third line: nominal values of the perihelion precessions due to the 
Sun's oblateness for J® = 2 x 10 -7 (Pireaux and Rozelot 2003; Pitjeva 2005a); the current level of uncertainty in it is about 
10% (Fienga et al. 2008). Fourth line: nominal values of the perihelion precessions due to the ring of the minor asteroids 
for m r i ng = 5 x 10~ 10 M© (Krasinsky et al. 2002); the uncertainty in it amounts to <5m r i ng = 1 x 10 -10 Mq (Krasinsky et 
al. 2002). Anyway, more recent estimates yield a smaller uncertainty: it is of the order of <5m r i ng = 0.1 — 0.3 x 10 -10 Mq 
(Fienga et al. 2009). Fifth line: nominal values of the perihelion precessions due to a massive ring modelling the action of 
the Classical Kuiper Belt Objects (CKBOs) for m CK BO s = 0.052 m e = 1.562 x 10 _7 M Q (Iorio 2007f). Actually, they may 
be smaller since direct estimates of the mass of the TNOs, modelled as a ring in the EPM2008 ephemerides, place an upper 
limit of 5.26 x 1O~ 8 M (Pitjeva 2010). 





Mercury 


Venus 


Earth 


Mars 


Azu 


-36 ± 50 [42] 


-4 + 5 [1] 


-2 + 4 [1] 


1 ± 5 [1] 




-20 


-3 


-1 


-0.3 


ZU ,0 
J 2 


+254 


+26 


+8 


+2 


ZU^ing 


+3 


+7 


+ 11 


+24 


MCKBOs 


+0.2 


+0.6 


+ 1 


+2 



(EPM2004) produced by Pitjeva (2005b). They are 
based on a data set of more than 317,000 observations 
(1913—2003) including radiometric measurements of 
planets and spacecraft, astrometric CCD observations 
of the outer planets and their satellites, and merid- 
ian and photographic observations. Such ephemerides 
were constructed by the simultaneous numerical inte- 
gration of the equations of motion for all planets, the 
Sun, the Moon, and the 301 largest asteroids. The ro- 
tations of the Earth and the Moon, the perturbations 
from the solar quadrupolar mass moment J® and that 
of the asteroid ring that lies in the ecliptic plane and 
which consists of the remaining smaller asteroids are 
included. With regard to the post-Newtonian dynam- 
ics, only the gravitoelectric, Schwarzschild-like terms 
of order 0(c -2 ), in the harmonic gauge (Newhall et al. 
1983), were included; the gravitomagnctic force of the 
Sun was not modelled. 

The EPM2004 ephemerides were used by Pitjeva 
(2005a) to phcnomenologically estimate corrections 
Aw to the known standard Newtonian/Einsteinian sec- 
ular precessions of the longitudes of perihelia of the in- 
ner planets as fitted parameters of a particular solution. 
Table 19 displays these corrections, which are obtained 
by comparing computer simulations based on the con- 
structed ephemerides with actual observations. In de- 
termining such extra-precessions the PPN parameters 
7 and (3 and the solar even zonal harmonic coefficient 
J® were not fitted; they were held fixed to their gen- 
eral relativistic and Newtonian values, i.e. 7 = ft = 1, 
J® = 2 x 10~ 7 . In the EPM2006 ephemerides (Pit- 
jeva 2008), used to accurately estimate the perihelion 



precession of Venus as well, were constructed by also 
modelling the actions of Eris and of the other 20 largest 
Trans- Neptunian Objects (TNOs); the database was 
enlarged by including, among other things, ranging 
data to the Magellan and Cassini spacecraft. 

Although the original purpose 31 of the estimation of 
the corrections Aw was not the measurement of the 
Lense-Thirring effect, the results of Table 19 can be 
used to take first steps towards an observational cor- 
roboration of the existence of the solar gravitomagnetic 
force. 

From Table 19 it turns out that the magnitude of 
the Lense-Thirring perihelion precessions of the inner 
planets just lies at the edge of the accuracy in determin- 
ing Aw. All the predicted Lense-Thirring precessions 
are compatible with the estimated corrections Aw. It 
should be noted that for Venus the 1 — a formal, statis- 
tical error in Aw is smaller than the gravitomagnctic 
effect. 



31 The goal of Pitjeva (2005a) was to make a test of the quality of 
the previously obtained general solution in which certain values 
of /3,7, J® were obtained. If the construction of the ephemerides 
was satisfactory, very small 'residual' effects due to such param- 
eters should have been found. She writes: "At present, as a 
test, we can determine [...] the corrections to the motions of 
the planetary perihelia, which allows us to judge whether the 
values of /3, 7 and J® used to construct the ephemerides are 
valid." . The smallness of the extra-perihelion precessions found 
in her particular test-solution is interpreted by Pitjeva as fol- 
lows: "Table 3 [of Pitjeva (2005a)] shows that the parameters 
7 = /3 = l,J = 2x 10 -7 used to construct the EPM2004 
ephemerides are in excellent agreement with the observations." 
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Table 20 Corrections Azi>, in 10~ 4 arcsec cy -1 (1 arcsec cy -1 = 10 mas yr -1 = 1.5 x 10~ 15 s -1 ), to the standard 
Newton/Einstein secular precessions of the perihelia of the inner planets estimated by Pitjeva (2010) with the EPM2008 
ephemerides including also radiometric data from Venus Express, Cassini and Magellan for Venus. The uncertainties are 
realistic. Source: Table 8 of Pitjeva (2010). Cfr. with Table 19. 



Mercury 


Venus 


Earth 


Mars 


Am -40 ± 50 


240 ± 330 


60 ±70 


-70 ± 70 



It must be considered that the inclusion of new data 
sets, modeling of further dynamical effects and new 
data processing techniques may have an impact on the 
estimation of the corrections Am. Thus, we are still in 
a preliminary, ongoing phase. As an example, in Table 
20, retrieved from Table 8 of Pitjeva (2010), we show 
new results for A-cc of the inner planets. 

Concerning the systematic alias due to the various 
competing dynamical effects listed in Table 19, the 
present-day level of their mismodelling would make 
them not particularly insidious (at ps 10% level of ac- 
curacy or even less, taking into account the recent im- 
provements in constraining the masses of the rings of 
the minor asteroids and of the TNOs), at least for Mer- 
cury and, especially, Venus, with the exception of the 
impact of the Sun's oblateness on Mercury. The latter 
could be removed anyway by suitably designing a lin- 
ear combination of the perihelia of Mercury and Venus, 
as done for the LAGEOS satellites. Moreover, an im- 
provement of the determination of J® by one order of 
magnitude with respect to the present-day level of un- 
certainty is one of the goals of the BcpiColombo mis- 
sion (Milani ct al. 2002; Ashby et al. 2007; Milani et al. 
2008). At present, it is known with an uncertainty of 
about 10% (Fienga et al. 2009). 

In principle, it would be possible to use also the 
nodes (Iorio 2005d), if only the corrections Afl to their 
standard precessions were available. Once other teams 
of astronomers have independently estimated their own 
corrections to the standard perihelion precessions (and, 
hopefully, to the nodal precessions as well) with differ- 
ent ephemerides, it will be possible to fruitfully repeat 
the present test. 

3.3 The interplanetary ranges 

An alternative approach to the perihelion precessions 
consists of looking at the Earth-planet ranges \p\, which 
are, among other things, directly observable quanti- 
ties 32 . We numerically investigated the solar Lense- 
Thirring effect on such an observable for the inner plan- 
ets (Iorio 2010d) by taking the difference between the 



perturbed range (|pp|) and the unperturbed, reference 
range (|ppj). The temporal interval of the numerical 
integrations has been taken equal to At = 2 yr because 
most of the present-day available time series of range 
residuals from several spacecraft approximately cover 
similar temporal extensions; moreover, also the typical 
operational time spans forecasted for future PLR tech- 
nique are similar. 

3.3.1 Earth-Mercury 

At present, the 1-way range residuals of Mercury from 
radar-ranging span 30 yr (1967-1997) and are at a few 
km-level (Figure B-2 of Folkner et al. (2008)); the same 
holds for the 1-way Mercury radar closure residuals cov- 
ering 8 yr (1989-1997, Figure B-3 a) of Folkner et al. 
(2008)). There are also a pair of Mariner 10 range resid- 
uals in the 70s at Mercury at 0.2 km level (Figure B-3 
b) of Folkner et al. (2008)). Ranging to BcpiColombo 
should be accurate to 4.5 — 10 cm (Ashby et al. 2007; 
Milani et al. 2008) over a few years. 

Figure 1 depicts the range perturbation due to 
the Sun's Lense-Thirring effect, neither considered so 
far in the dynamical force models of the planetary 
ephemerides nor in the BepiColombo analyses so far 
performed. The peak-to-peak amplitude of the Lense- 
Thirring signal is up to 17.5 m over 2 yr, which, if 
on the one hand is unmeasurable from currently avail- 
able radar-ranging to Mercury, on the other hand cor- 
responds to a potential relative accuracy in measuring 
it with BepiColombo of 2 — 5 x 10~ 3 ; this clearly shows 
that the solar gravitomagnetic field should be taken into 
account in future analyses and data processing. Other- 
wise, it would alias the recovery of other effects. 

Figure 2 shows the nominal signature of the Sun's 
quadrupolar mass moment on the Mercury range for 
jf = 2x 10~ 7 . Its action has been modeled as (Vrbik 
2005) 



3J°i?|GM G 
2r 4 



5 ( f-k 



(f-k) 



(49) 



32 Recall that the perihelia, like all the other Keplerian orbital 
elements, cannot be directly measured. 



where Rq is the Sun's mean equatorial radius and k is 
the unit vector of the z axis directed along the body's 
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EMB-Mercury range: SSB numerical calculation 
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Fig. 1 Difference A\p\ = \pp\ — \p\\ in the numerically in- 
tegrated EMB-Mercury ranges with (pp) and without (pp.) 
the perturbation due to the Sun's Lense-Thirring field over 
At = 2 yr. The same initial conditions (J2000) have been 
used for both integrations. The state vectors at the refer- 
ence epoch have been retrieved from the NASA JPL Hori- 
zons system. The integrations have been performed in the 
ICRF/J2000.0 reference frame, with the mean equinox of 
the reference epoch and the reference {xy} plane rotated 
from the mean ecliptic of the epoch to the Sun's equator, 
centered at the Solar System Barycenter (SSB). 

rotation axis. The signal in Figure 2 has a maximum 
span of 300 m, corresponding to an accuracy measure- 
ment of 3 x 10~ 4 . A measure of the solar jf accurate 
to 10 -2 is one of the goals of BcpiColombo (Milani et 
al. 2002); knowing precisely jf would yield important 
insights on the internal rotation of the Sun. 

Concerning the impact of neglecting the gravitomag- 
netic field of the Sun in future data analyses, it may 
affect the determination of jf at the 12% level. On 
the other hand, in order to allow for a determination 
of the Lense-Thirring effect, the Sun's quadrupolc mass 
moment should be known with an accuracy better than 
at present by at least one order of magnitude; this is 
just one of the goals of BepiColombo. Anyway, also the 
time signatures of the two signals would play a role. 

3.3.2 Earth-Venus 

Although at present, contrary to Mercury and Mars, 
no tests of interplanetary ranging to Venus have been 
practically performed, we include this case not only 
for completeness but also because simulations of inter- 
planetary transponder and laser communications exper- 
iments via dual station ranging to SLR satellites cover- 
ing also Venus have been implemented (Degnan 2006, 
2008). Currently available radar-ranging normal points 
to Venus cover about 33 yr, from 1962 to 1995. The 
range residuals are depicted in Figure B-6 of Folkner et 
al. (2008); after having been as large as 15 km in the 
first 10 yr, they have dropped below 5 km since. Figure 



EMB-Mercury range: SSB numerical calculation 
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Fig. 2 Difference A\p\ = \pp\ — |pp, | in the numerically in- 
tegrated EMB-Mercury ranges with (pp) and without (pn) 
the nominal perturbation due to the Sun's quadrupole mass 
moment J® = 2.0 x 10~ 7 over At — 2 yr. The same ini- 
tial conditions (J2000) have been used for both integrations. 
The state vectors at the reference epoch have been retrieved 
from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, 
with the mean equinox of the reference epoch and the refer- 
ence {xy} plane rotated from the mean ecliptic of the epoch 
to the Sun's equator, centered at the Solar System Barycen- 
ter (SSB). 

B-4 of Folkner et al. (2008) shows the range residuals 
to Venus Express at Venus from 2006 to 2008; the are 
below the 10 m level. 

Figure 3 shows the Lense-Thirring perturbation of 
the Venus range, integrated in a frame aligned with 
the Sun's equator. The peak-to-peak amplitude is 2 m, 
which would be measurable with a future accurate cm- 
level ranging device with a relative accuracy of 2 — 5 x 
10~ 2 . The Lense-Thirring signature is still too small 
to be detected nowadays with the current spacecraft 
ranging. 

The range perturbation due to the Sun's oblateness 
is depicted in Figure 4 for the nominal value J® — 
2 x 10~ 7 . Also in this case a barycentric frame rotated 
to the Sun's equator has been adopted. The nominal 
maximum shift is 40 m, so that a measure accurate to 
2.5 x 10~ 3 would be possible with a future 10 cm-level 
ranging technique. 

If not modeled, the Sun's gravitomagnetic field 
would impact a determination of J® at the 5% level. 
Conversely, the present-day 10% uncertainty in the 
Sun's oblateness would yield a mismodeled signal two 
times larger than the gravitomagnetic one. However, 
their temporal signatures are different and would there- 
fore allow to seperate these two effects. 
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Earth- Venus range: SSB numerical calculation 
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Fig. 3 Difference A|p] = pp| — |pr| in the numerically 
integrated EMB- Venus ranges with (pp) and without (pn) 
the perturbation due to the Sun's Lense-Thirring field over 
At = 2 yr. The same initial conditions (J2000) have been 
used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL 
Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the mean equinox 
of the reference epoch and the reference {xy} plane rotated 
from the mean ecliptic of the epoch to the Sun's equator, 
centered at the Solar System Barycenter (SSB). 

Earth-Venus range: SSB numerical calculation 
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Fig. 4 Difference A\p\ = |pp — |pR in the numerically 
integrated EMB- Venus ranges with (pp) and without (pr) 
the nominal perturbation due to the Sun's quadrupole mass 
moment J® = 2.0 x 10 -7 over At = 2 yr. The same initial 
conditions (J2000) have been used for both the integrations. 
The state vectors at the reference epoch have been retrieved 
from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, 
with the mean equinox of the reference epoch and the refer- 
ence {xy} plane rotated from the mean ecliptic of the epoch 
to the Sun's equator, centered at the Solar System Barycen- 
ter (SSB). 

3.3.3 Earth-Mars 

For Mars we have at our disposal long time series of 
range residuals accurate to about 1 — 10 m-level thanks 
to several spacecraft (Viking, Mars Pathfinder, Mars 



Earth-Mars range: SSB numerical calculation 
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Fig. 5 Difference A\p\ = \pp\ — \pn\ in the numerically 
integrated EMB-Mars ranges with (pp) and without (pk) 
the perturbation due to the Sun's Lense-Thirring field over 
At = 2 yr. The same initial conditions (J2000) have been 
used for both integrations. The state vectors at the refer- 
ence epoch have been retrieved from the NASA JPL Hori- 
zons system. The integrations have been performed in the 
ICRF/J2000.0 reference frame, with the mean equinox of 
the reference epoch and the reference {xy} plane rotated 
from the mean ecliptic of the epoch to the Sun's equator, 
centered at the Solar System Barycenter (SSB). 

Global Surveyor, Mars Odyssey, Mars Reconnaissance 
Orbitcr, Mars Express) which have orbited or are still 
orbiting the red planet. Figure B-10 of (Folkner et al. 
2008) depicts the 1-way range residuals of the Viking 
Lander at Mars spanning from 1976 to 1982; they 
are approximately at the 20 m level. Figure B-ll of 
(Folkner et al. 2008) shows the 1-way range residuals 
of several post- Viking spacecraft; they generally cover 
a few years and are accurate to 5 — 10 m. 

The Lense-Thirring range perturbation, computed 
in a frame aligned with the Sun's equator, is shown in 
Figure 5. Its peak-to-peak amplitude is about 1.5 m, 
not too far from the present-day range accuracy; thus, 
its existence as predicted by general relativity is not 
in contrast with the range residuals currently available. 
It could be measured with a future cm-level ranging 
system at a 3 — 6% level. 

Figure 6 illustrates the nominal signal due to the 
Sun's quadrupole mass moment for jf — 2 x 10~ 7 
computed in a frame aligned with the Sun's equator. 
Its peak-to-peak amplitude amounts to 25 m; thus, its 
effect would be well measurable at a 2 — 4 x 10 -3 level 
by means of a new ranging facility with an accuracy 
of the order of cm. Since J® is nowadays accurate to 
10 _1 , the corresponding mismodeled signature would 
be as large as about 2.5 m. If not properly modeled, 
the Lense-Thirring effect would bias the J® signal at 
a 6% level. Conversely, considering J® as a potential 
source of systematic bias for the recovery of the gravito- 
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Earth-Mars range: SSB numerical calculation 
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Fig. 6 Difference A\p\ = \pp\ — \pr\ in the numerically 
integrated EMB-Mars ranges with (pp) and without (pn) 
the nominal perturbation due to the Sun's quadrupole mass 
moment J® = 2.0 x 10~ 7 over At — 2 yr. The same ini- 
tial conditions (J2000) have been used for both integrations. 
The state vectors at the reference epoch have been retrieved 
from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, 
with the mean equinox of the reference epoch and the refer- 
ence {xy} plane rotated from the mean ecliptic of the epoch 
to the Sun's equator, centered at the Solar System Barycen- 
ter (SSB). 

magnetic effect, the mismodclcd signature of the Sun's 
quadmpolar mass moment would be 1.6 times larger 
than it. An improvement in its knowledge by one or- 
der of magnitude, as expected from, e.g., BepiColombo, 
would push its bias on the Lense-Thirring signal at 16%. 
However, it must be noted that their temporal evolu- 
tions arc different. 



4 The gravitomagnetic field of Mars: 
perspectives for its detection 

5 General overview 

Since the angular momentum of Mars can be evaluated 
to be 

S M = (1-92 ± 0.01) x 10 32 kg m 2 s" 1 (50) 

from the latest spacecraft-based determinations of the 
areophysical parameters (Konopliv ct al. 2006), the cor- 
responding gravitomagnetic length reads 



This has to be compared with the present-day accuracy 
in determining the orbit of a spacecraft like, e.g., Mars 
Global Surveyor (MGS) which is about 0.15 m in the ra- 
dial direction (Konopliv et al. 2006) and not affected by 



the gravitomagnetic force itself. Thus, it makes sense, 
in principle, to investigate the possibility of measuring 
the Lense-Thirring effect in the gravitational field of 
Mars as well. 

In fact, Iorio (2006c, 2007g) proposed an interpreta- 
tion of the time series of the RMS orbit overlap differ- 
ences (Konopliv et al. 2006) of the out-of-plane part v of 
the orbit of MGS over a time span AP of about 5 years 
(14 November 1999-14 January 2005 in Iorio (2007g)) 
in terms of the Lense-Thirring effect. It turned out 
that the average of such a time series over AP, normal- 
ized to the predicted Lense-Thirring out-of-planc mean 
shift over the same time span, is /i = 1.0018 ± 0.0053. 
The interpretation by Iorio (2006c, 2007g) has recently 
been questioned by Krogh (2007); a reply has been set 
in Iorio (2010e). Basically, various linear fits to dif- 
ferent data sets including, among others, the full time 
series of the entire MGS data (4 February 1999-14 Jan- 
uary 2005) were made as well; the predictions of general 
relativity turn out to be always confirmed. The analyt- 
ical calculation of the competing aliasing effects due to 
both the gravitational and non-gravitational perturba- 
tions, which affect the in-plane orbital components of 
MGS, do not show up in the real data. Moreover, the 
non-conservative forces, whose steadily refined model- 
ing mainly improved the in-plane orbital components 
of MGS, not the normal one, exhibit high-frequency, 
non-cumulative in time variations (Forbes et al. 2006). 

In view of the never fading interest for planetologi- 
cal missions to the red planet, the preliminary design 
of a spacecraft-based dedicated mission to Mars has 
been investigated by Iorio (2009b). Here we will deal 
with such an issue in some detail. In particular, we 
will concentrate on the the multipolar expansion of the 
martian gravitational potential in order to see what are 
the critical issues in view of the present-day knowledge 
of the Martian space environment. We will not dis- 
cuss here the perturbations of non-gravitational origin 
which depend on the shape, the instrumentation and 
the orbital maneuvers of such probes. They would be 
strongly related to possible other tasks, more consistent 
with planetology, which could be fruitfully assigned to 
such a mission in order to enhance the possibility that 
it may become something more than a mere, although- 
hopcfully-interesting, speculation; in this respect the 
medium-long term ambitious programs of NASA to 
Mars may turn out to be useful also for the purpose 
discussed here. 

5.1 The use of one nearly polar spacecraft 

Let us start with a polar orbital geometry and exam- 
ine the systematic error 5ji induced by the uncertainty 
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Fig. 7 Systematic percent errors Sfij e per degree £ due 
to the mismodelling in the even zonal harmonics SJi,£ = 
2, 4, 6... for a = 3796 km, e = 0.01, 89 deg < I < 91 deg ac- 
cording to the calibrated covariance sigmas of the MGS95J 
global gravity solution up to £ — 20. The Lense-Thirring 
effect amounts to 34 mas yr _1 corresponding to a shift in 
the out-of-plane direction of 0.62 m yr _ . 



a = 9500 km, e = 0.01 
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Fig. 8 Systematic percent errors S/j,j e per degree £ due 
to the mismodelling in the even zonal harmonics 5Ji,£ = 
2,4,6... for a = 9500 km, e = 0.01, 89 deg < / < 91 
deg according to the calibrated covariance sigmas of the 
MGS95J gravity model up to £ = 20. The Lense-Thirring 
effect amounts to 2 mas yr _1 corresponding to a shift in the 
out-of-plane direction of 0.10 m yr~ . 



in the even zonals 5Ji on the node VL by assuming 
for them the calibrated covariance sigmas of the latest 
Mars gravity model MGS95J (Konopliv et al. 2006). 
We will evaluate it as 

8nj t <\tl.t\8Jt,l = 2,A,--. (52) 

Because of the presence of cos / in all the coefficients 
of the oblateness-induced node precessions, we will 
concentrate on nearly polar (/ sa 90 deg) orbital config- 
urations to minimize such a corrupting effect. It turns 
out that altitudes of a few hundreds of km typical of the 
majority of the currently ongoing martian missions are 
definitely not suited for our scope, as shown by Figure 
7. It must also be pointed out that likely more even 
zonals of degree higher than I = 20 would come into 
play for such low orbital altitudes, as in the case of the 
forthcoming terrestrial LARES satellite. Note that in- 
serting a probe into an areocentric orbit is not an easy 
task, so that we decided to allow for a departure of up 
to 1 deg from the ideal polar orbital configuration to 
account for unavoidable orbital injection errors. Also 
typical mission requirements pull the inclination some 
degrees apart from 90 deg: for Mars Global Surveyor 
(MGS) I = 92.86 deg. The systematic bias per degree 
increases for higher degrees I = 2, 4, 6, . . . and it turns 
out that only a very narrow range for /, i.e. AI ss 10~ 3 
deg, unlikely to obtain, might push the systematic er- 
rors below the 1% level. 

By keeping a near polar geometry, much larger values 
of the scmimajor axis, comparable with that of Phobos 
(a = 9380 km) , one of the two natural satellites of Mars, 



yield reasonable results, as shown in Figure 8. Indeed, 
now the noise decreases with high degree terms and the 
most serious contribution is due to the first even zonal 
I = 2 yielding a bias up to about 11% for AI = 1 deg; J4 
affects the the Lense-Thirring precessions at most at the 
3% level for AI = 1 deg. Note also that a much larger 
departure from 90 deg (AI rj 0.1 deg) would allow a 
further reduction of the noise below 1%. It must be 
noted that for such a high-altitude orbital configuration 
the gravitomagnetic shift in the out-of-plane direction 
would amount to 0.10 m yr ; it is certainly a small 
figure, but, perhaps, not too small if one considers that 
the average shift in the out-of-plane direction of the 
much lower Mars Global Surveyor is 1.6 m after about 
5 yr. It poses undoubtedly challenging requirements 
in terms of sensitivity and overall orbit determination 
accuracy, but future improvements may allow to detect 
such a small displacement. 

Finally, let us remark that, in principle, also the tem- 
poral variations of J2 should be accounted for; however, 
since at present no secular trends have been detected, 
such changes, mainly seasonal, annual and semi-annual 
(Konopliv et al. 2006) would not seriously impact our 
measurement. 

Another martian parameter which must be taken 
into account is the equatorial radius R along with its 
uncertainty. Since 

Cl. e oc R l , (53) 
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a = 9500 km, e = 0.01, SR = 0.04km 
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Fig. 9 Systematic percent errors Sfij e per degree £ due 
to the uncertainty in Mars' radius R, assumed to be SR = 
6080 m, for a = 9500 km, e = 0.01, 89 deg < i < 91 deg. 
The errors for I = 4, 6 are as large as 6000% and 500%, 
respectively. 



the systematic error per degree induced by SR can be 
written as 
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(54) 



If we assume conservatively for it SR = 6080 m, i.e. 
the difference between the reference value of MGS95J 
(Konopliv et al. 2006) and the one in Yoder (1995), the 
result is depicted in Figure 9 for a = 9500 km. If, in- 
stead, one takes 8R = 0.04 km the errors per degree 
are as in Figure 10 for a = 9500 km. Obviously, R 
is the most serious source of systematic error, even for 
high altitudes. Moreover, it might be unrealistic to ex- 
pect improvements in the determination of R by future 
areocentric missions able to push SR below terrestrial 
values, i.e. 1 — 0.1 m. Thus, R will likely remain an 
insurmountable obstacle if only one probe is to be used. 

Concerning GM, it turns out that the impact of its 
mismodelling is of no concern being < 1%. 

In summary, in this Section we investigated a nearly 
polar, high-altitude Phobos-like orbit and noted that 
the Lense-Thirring effect amounts to a shift of 0.10 m 
yr _1 in the out-of-plane direction. Concerning the sys- 
tematic errors, the most crucial source of aliasing is the 
radius R (Sfi w 10000% for SR = 0.04 km) and, to a 
lesser extent, the first even zonal harmonic J 2 of the 
areopotential (S[i < 11%). To reduce them the probe 
should be inserted into an orbit with an inclination close 
to 90 deg within 10~ 3 deg or less. In conjunction with 
such a very tight constraint, it must also be hoped that 
future missions to Mars will improve our knowledge of 
the fundamental parameters of the red planet to a suffi- 
cient extent to allow for larger departures from the ideal 



Fig. 10 Systematic percent errors 5^j e per degree £ due to 
the uncertainty in Mars' radius R, assumed to be 5R = 0.04 
km, for a = 9500 km, e = 0.01, 89 deg < i < 91 deg. The 
errors for £ = 4, 6 are as large as 40% and 5%, respectively. 



polar geometry: after all, the gravity solution MGS95J 
represents an improvement of one order of magnitude 
with respect to the previous MGS75D model (Yuan et 
al. 2001). However, this may be valid for the spherical 
harmonic coefficients of the areopotential, not for the 
radius of Mars for which an uncertainty of the order 
of one meter would not be enough. We must conclude 
that the use of only one probe for measuring its nodal 
Lense-Thirring precession is unfeasible. 

5.2 Including Phobos and Deimos 

A possible solution to reduce the systematic errors 
would be to consider a linear combination of the nodes 
of the proposed probe and of the two natural satellites 
of Mars, i.e. Phobos (a = 9380 km, I = 1.075 deg, 
e = 0.0151) and Deimos (a = 23460 km, I = 1.793 deg, 
e = 0.0002), suitably designed to cancel out the impact 
of just SJ2 and SR/R, according to an approach sug- 
gested for the first time in the context of the terrestrial 
LAGEOS-LAGEOS II test (Ciufolini 1996). A possible 
combination is 
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with Sti denoting an Observed-minus-Calculated (O — 
C) quantity accounting for unmodcllcd/mismodelled 
features of motion and extracted from the data pro- 
cessing with full dynamical force models. By purposely 
leaving the gravitomagnctic force unmodelled, it can be 
written as 



Sfl = ^A LT + n. 2 SJ 2 + SI.r 



SR 
~R 



(56) 
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here /i = 1 in general relativity and the coefficient £1,r 
is defined as 



e=2 



dR 



(57) 



and A represents all other remaining mismodcllcd/un- 
modcllcd effects acting on the node (e.g. 8J4, 5Jq, . . . ; 
SGM). By solving for /it one gets 



SPhobosADc 



Oil 



a 2 



(58) 



Note that «i, contrary to a.2, is not defined for I = 90 
deg because f_^ robo , which is proportional to cos J, van- 
ishes for all £ at I = 90 deg. Now, the systematic bias 
of J4, J6, Js r ■ • on the combination of cq. (55) turns 
out to be Sfi < 6% for 89 deg < I < 89.9 deg where 
the major contribution is due to J 4. It is an acceptable 
result, especially in view of the fact that the admissi- 
ble range for the inclination amounts now to about 1 
deg and that eq. (55) is, by construction, immune to 
the uncertainty in the martian radius; indeed, although 
it is likely that many physical properties and parame- 
ters like, e.g., the even zonals, of the red planets will 
be determined with increasing accuracy by the many 
ongoing and planned missions, it is unlikely that the 
radius will be known to a sufficient accuracy to change 
at an acceptable level the bias induced by it (see Figure 
9 and Figure 10). 

The combination of eq. (55) is also able to remove a 
large part of the effect of the mismodclling in GM 



v 7 e=2 
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(59) 



Indeed, the main contribution to eq. (59) is due to J 2 , 
which is canceled out by eq. (55) with the coefficients 
of eq. (58). The other terms in cq. (59), not canceled 
by the combination of eq. (55), yield negligible errors 
well below 1%. 

However, it must be stressed that the node of Phobos 
undergoes secular precessions due to other perturba- 
tions of gravitational origin (Lainey et al. 2007) (e.g., 
the non-sphericity of Phobos itself, the martian tidal 
bulge and nutation) which have not been considered 
here and that would affect the combination of eq. (55). 
It turns out that the most relevant one is that due to 
the spherical harmonic coefficients C20 and C22 of Pho- 
bos itself (Bordcrics and Yoder 1990) whose induced 
secular precession amounts nominally to 200 km over 
3 yr (Lainey et al. 2007). The nutation perturbation 



nominally amounts to 0.3 km over 3 yr (Lainey et al. 
2007). The tidally induced precession, parameterized 
in terms of the martian Love number k 2 would amount 
nominally to 0.06 km after 3 yr (Lainey et al. 2007), 
but being k 2 = 0.152 ±0.009 (Konopliv et al. 2006), its 
modelling would left a w 1 m yr -1 mismodelled trend. 
The Phobos gravitomagnetic shift is 0.1 m yr -1 . 

In terms of the detectability of the Lense-Thirring 
signal with the combination of eq. (55), undoubtedly a 
major drawback of the strategy of including Phobos and 
Deimos is the poor accuracy with which their orbits can 
be reconstructed with respect to their Lense-Thirring 
signal. Indeed, while their gravitomagnetic shifts are of 
the order of 1 — 10 cm yr -1 , the latest NASA martian 
ephemerides 33 for them yield an accuracy of 1 — 10 km 
in the radial, transverse and out-of-plane orbital com- 
ponents (Jacobson and Rush 2006). Thus, according 
to the present-day level of accuracy, including the nat- 
ural satellites of Mars in the combination of eq. (55) 
would introduce a noise which would overwhelm the 
relativistic trend of interest. This situation may be- 
come more favorable in the near future in view of the 
planned Phobos-Grunt space mission 34 (Marov et al. 
2004) which should allow, among other things, to im- 
prove the orbit determination of Phobos as well. 

5.3 Two dedicated probes 

In principle, the linear combination approach could be 
followed by using the nodes of two other spacecraft 35 , 
although sending to Mars three new probes would in- 
crease the costs and the difficulties of such a demanding 
mission. 

It is interesting to consider a scenario involving 
only two probes, named PI and P2, in a nearly po- 
lar counter-orbiting configuration, proposed for the first 
time by Van Patten and Everitt (1976a) in the frame- 
work of the attempts to design a suitable terrestrial 
mission, and later generalized by Ciufolini (1986) for 
other inclinations. By assuming for, say, PI ai = 9500 
km, 7i = 89 deg, e\ = 0.01, it turns out that the sum 
of their nodes would be an observable relatively insen- 
sitive to departures from the ideal configuration for P2, 
i.e. a 2 = 9500 km, I 2 = 91 deg, e 2 = 0.01, at least 
in regard to the mismodelling in the even zonals. For 
a 2 = 9600 km, 90.05 deg < I 2 < 92 deg, e 2 = 0.03 
the bias due to J 2 is up to 7%, while the other even 

33 See on the WEB http://ssd.jpl.nasa.gov/7sat_ephem. 

34 See on the WEB http://www.esa.int/esaMI/ESA_Permanent_ 
Mission jn_Russia/SEMIJFW4QWD_0.html. 

35 For a polar geometry the Lense-Thirring precession of the peri- 
centre u of a spacecraft vanishes being proportional to cos /. 
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5.4 Summary 




Fig. 11 Systematic errors 8fi,j e per degree due to the mis- 
modelling SJ2, SJa, 8J§, ... (MGS95J model) in the uncan- 
celled even zonal harmonics in the combination of the nodes 
of probe PI (ai = 9500 km, h = 89 deg, ei = 0.01) and 
probe P2 (a 2 = 9600 km, 60 deg < h < 89 deg, e 2 = 0.03) 
with eq. (61) which cancels out entirely the bias due to SR. 



zonals would have an impact of the order of 1% or less. 
Unfortunately, the uncertainty in the martian radius 
makes such constraints much more stringent. Indeed, 
S^r ps 10 5 % for SR = 6080 m and 6fi R ps 4000% for 
SR = 0.04 km. Also in this case, the improvement in 
R expected from future missions to the red planet may 
be not sufficient. 

An approach that could be successfully followed 
within the framework of the linear combination strat- 
egy with two probes PI and P2, especially in view of 
likely future improvements in the even zonal harmonics 
of the areopotential, consists in designing a combina- 
tion which, by construction, entirely cancels out the 
bias due to the uncertainty in R 



SR 



(60) 



being, instead, affected by SJ2, 6J4, SJq. It is analogous 
to the three- nodes combination of eq. (32), with 



(61) 



Indeed, as can be seen from Figure 11, Sixj t would 
be rather small (< 1%), according to the present-day 
MGS95J model, provided that inclinations I <C 90 deg 
are adopted for the second probe. It turns out that 
the most relevant even zonals would be J4 and Jq with 
a bias of the order of ss 1 — 0.5%. With eq. (61) an 
inclination range as large as ten degrees would be ad- 
missible: this is really a great advantage with respect 
to the other cases examined so far. 



The main source of systematic bias is the mismodclling 
in some of the parameters entering the multipolar ex- 
pansion of the Newtonian part of the areopotential, es- 
pecially the Mars' equatorial radius R and, to a lesser 
but non- negligible extent, the even zonal harmonics 
Je, I = 2, 4, 6, Since the resulting aliasing node pre- 
cessions are proportional to R e a~( 3 / 2+e ^ cos/, a high- 
altitude (a ps 9500 km), polar (/ = 90 deg) probe would 
be an optimal solution, but the estimation of the un- 
avoidable orbit injection errors showed that such an op- 
tion is unfeasible because it would require too stringent 
constraints on the departures from the ideal case / = 90 
deg. In principle it would be possible to suitably com- 
bine the nodes of such a probe with those of the natural 
satellites of Mars, Phobos and Deimos. However, the 
Lense-Thirring shifts of these bodies (1 — 10 cm yr _1 ) 
are orders of magnitude too small with respect to the 
present accuracy in reconstructing their orbits (1 — 10 
km). 

A viable option consists in using two probes PI 
and P2 at high-altitudes (01 = 9500 km, a 2 = 9600 
km) and different inclinations (Ji ps 90 deg, 60 deg 
< I2 < 89 deg), and combining their nodes so as to 
entirely cancel out the bias due to R: the resulting bias 
due to SJ2,SJi,SJe,... would be < 1%, according to 
the present-day MGS95J gravity model, over a range 
of values for the inclination as large as ten degrees. 

A major challenge would certainly be to reach a sat- 
isfying accuracy in reconstructing the orbits of such 
probes whose Lense-Thirring out-of-plane shifts would 
amount to about 10 cm yr" 1 ; for example, the space- 
craft should remain operative for many years, without 
any failure in the communication with the Earth, and 
it is likely that one or more landers would be required 
as well. 

Finally, we emphasize that a different type of anal- 
ysis should be done to complete the studies presented 
here, based on direct observable quantities like range 
and range-rates. 
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6 The gravitomagnetic field of Jupiter: 

perspectives in measuring the jovian angular 
momentum via the Lense-Thirring effect 

6.1 General overview 

After the original proposal by Lense and Thirring 
(1918) of looking at the Galilean satellites 36 of Jupiter 
as natural probes for measuring the gravitomagnetic 
precessions of their longitudes of perijove w. Bogorod- 
skii (1959) considered the centennial shift of the ar- 
gument of perijove u> of the inner moon Amalthea 
(a = 181400 km, e = 0.0032) amounting to 112 arc- 
sec, while Soffcl (1989) looked at centennial node shift 
of Io of about 9 arcsec. No evaluations of the system- 
atic bias posed by several competing dynamical effects 
of the Jovian space environment have been examined 
in such earlier studies. In more recent times, Iorio 
and Lainey (2005) revisited the Jovian system of the 
Galilean moons in view of the recent advancements in 
their orbit determination. It turned out that most of 
their gravitomagnetic signals would be absorbed in the 
estimation procedure of the satellites' initial state vec- 
tors; the remaining signatures would amount to some 
tens of meters, while the current best observations 
have an accuracy of a few tens of kilometers 37 . Iorio 
and Lainey (2005) concluded that a spacecraft orbiting 
Jupiter would be required to reveal its gravitomagnetic 
field; as we will see below, such an opportunity may be 
offered in the near future (Iorio 2010f). 

In an attempt to measure gravitomagnctism by 
means of the deflection of electromagnetic waves by 
Jupiter due to its orbital motion 38 , a dedicated VLBI- 
bascd radio-interferometric experiment has been per- 
formed (Kopeikin & Fomalont 2006; Fomalont and 
Kopcikin 2008). With regard to other suggested non- 
gravitomagnetic tests of general relativity in the jovian 
gravitational field, Hiscock and Lindblom (1979) pro- 
posed to measure the much larger gravitoelectric Ein- 
stein pericenter precessions of the natural satellites of 
Jupiter and Saturn. There exist also plans for perform- 
ing a test of the light bending due to Jupiter's monopole 
and quadrupole mass moments with the forthcoming 
astrometric mission GAIA (Crosta and Mignard 2006). 
Finally, also the proposal by Haas and Ross (1975) to 



36 Their semimajor axes range from 421800 km (Io) to 1882700 
km (Callisto), the eccentricities being of the order of 10 — 3 . See 
on the WEB http://ssd.jpl.nasa.gov/7sat_elem. 

37 See on the WEB http://ssd.jpl.nasa.gOv/7sat_ephem#ref3. 

38 In this case, the mass currents inducing a gravitomagnetic ac- 
tion are not those related to the Jupiter's proper rotation (in- 
trinsic gravitomagnctism), but arc due to its translational orbital 
motion (extrinsic gravitomagnctism). 



Table 21 Planetocentric nominal orbital parameters of 
Juno, a, e, / are the semi-major axis (in jovian radii R = 
71492 km) , the eccentricity and the inclination (in deg) with 
respect to Jupiter's equator, respectively. Qlt and cjlt are 
Juno's Lense-Thirring node and perijove precessions (in mas 
yr" 1 ) for S = 6.9 x 10 38 kg m 2 s _1 (Soffel et al. 2003). P 
is the orbital period in days. T is the mission duration in 
years. 
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measure Schiff's gyroscope precession in the gravita- 
tional field of Jupiter should be mentioned. 

Another opportunity to fruitfully use the Jovian 
system may open after the approval of the Juno mis- 
sion 39 by NASA (Matousek 2007) to Jupiter. Juno is 
a spinning, solar powered spacecraft to be placed in 
a highly eccentric polar orbit around Jupiter (see Ta- 
ble 21 for its relevant orbital parameters) specifically 
designed to avoid its highest radiation regions. Un- 
derstanding the formation, evolution and structure of 
Jupiter is the primary science goal of Juno. It will 
carry onboard a dual frequency gravity/radio science 
system, a six wavelength microwave radiometer for at- 
mospheric sounding and composition, a dual-technique 
magnetometer, plasma detectors, energetic particle de- 
tectors, a radio/plasma wave experiment, and an ul- 
traviolet imager/spectrometer. The nominal mission's 
lifetime is 1 year. Juno is aimed, among other things, 
at accurately mapping the gravitational field of Jupiter 
(Anderson 1976) with unprecedented accuracy (Ander- 
son et al. 2004) by exploiting the slow apsidal precession 
of its 11-day orbit. 

Since the Lense-Thirring precession is due to the 
proper angular momentum S of the orbited central 
body, one may also take its existence as granted and 
use it as a direct, dynamical measurement of Jupiter's 
angular momentum through the Lense-Thirring effect; 
this would yield further, important information con- 
cerning the interior of Jupiter. Indeed, the moment 
of inertia ratio C /MR 2 entering S is a measure of the 
concentration of mass towards the center of the planet 
(Irwin 2003). Such a figure, together with the mea- 
sured values of the zonal 40 coefficients of the gravity 
field accounting for its deviations from spherical sym- 
metry may be fitted with internal models of the vari- 
ation of the density, pressure, temperature and com- 
position with depth (Irwin 2003; Guillot 2005; Hori et 
al. 2008). Moreover, a dynamical, model-independent 



See on the WEB http://juno.wisc.edu/index.html 
They preserve the axial symmetry. 
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determination of S would be important also for a bet- 
ter knowledge of the history and formation of Jupiter 
(Machida et al. 2008). 

6.2 Sensitivity analysis 

Jupiter's proper angular momentum amounts to (Soffcl 
et al. 2003) 

5sj 6.9 x 10 38 kg m 2 s" 1 (62) 

and the corresponding Lense-Thirring precessions of 
Juno's node and perijove are listed in Table 21 . It 
turns out that only the node experiences a GM preces- 
sion corresponding to a shift Av of the cross-track com- 
ponent of the planetocentric position (Christodoulidis 
et al. 1988) 



Ai/ LT = a\j 1 + — sin IAQ LT = 572 m (T = 1 yr) (63) 

over the entire duration of the mission. A total accu- 
racy of the order of 1-10 m with respect to the km-level 
of the past Jupiter missions in reconstructing Juno's 
orbit in a planetocentric frame does not seem an un- 
realistic target, although much work is clearly required 
in order to have a more firm answer. Note that a 1-10 
m accuracy implies a 0.2 — 2% error in measuring the 
gravitomagnctic shift 

To consider the possible detection of the Lense- 
Thirring effect by means of the Juno mission is also sug- 
gested by a different approach with respect to the cumu- 
lative measurement over the full mission duration previ- 
ously outlined. Indeed, a gravity-science pass for Juno 
is defined by a continuous, coherent Doppler range-rate 
measurement plus and minus three hours of closest ap- 
proach; in practice, most of the Lense-Thirring preces- 
sion takes place just during such a six-hours pass, a 
near optimum condition. Another crucial factor is the 
orientation of the Earth with respect to Juno's orbit: 
our planet will be aligned 67 deg from the probe's or- 
bital plane at approximately two degrees south latitude 
on the jovian equator. Preliminary numerical simula- 
tions of Juno's Lense-Thirring Doppler range-rate sig- 
nal show that such an orbital geometry represents a 
perfect compromise for measuring both Jupiter's even 
zonal harmonics and the gravitomagnctic signal itself. 
Indeed, it turns out that the maximum Lense-Thirring 
Doppler signal over a single six-hour gravity pass is of 
the order of hundred /im s _1 , while the limit of accu- 
racy for Juno's Ka-band Doppler system is about one 
/urn s _1 over such a pass. Thus, even by taking 25 
repeated passes out of a total of approximately 33, it 



would be possible to reach a measurement precision be- 
low the percent level. 

It should be noted that there are estimates for S in 
the literature in favor of smaller values than given in 
eq. (62) by a factor 1.5 — 1.6; for example, Machida et 
al. (2008) report 

S = 4.14 x 10 38 kg m 2 s _1 ; (64) 

the ratio of of the values of eq. (62) to cq. (64) 
is 1.6, i.e. close to 1.5 coming from the ratio of 
C/MR 2 = 2/5 = 0.4, valid for a homogenous sphere, to 
CI MR 2 = 0.264 by Irwin (2003) who assumes a con- 
centration of mass towards Jupiter's center. Here we 
consider only the systematic uncertainty induced by the 
imperfect knowledge of the Newtonian part of Jupiter's 
gravitational field and use the value of eq. (62) for S. 

Also in this case, a major source of systematic uncer- 
tainty is due to the deviation of Jupiter's gravitational 
field from spherical symmetry (Anderson 1976). 

6.2.1 Analytical calculations 

As seen in Section 2.1, the zonal (m = 0) harmonic coef- 
ficients Ji of the multipolar expansion of the Newtonian 
part of the planet's gravitational potential give rise to 
long-period, (i.e. averaged over one orbital revolution), 
orbital effects on the longitude of the ascending node 
O, the argument of pericentre ui and the mean anomaly 
Ai of the form 

(*) =^2^fuJt, * = Sl,u,M, (65) 
e=2 

where are coefficients depending on the planet's 
GM and equatorial radius R, and on the spacecraft's 
inclination / and eccentricity e through the inclination 
Fimp(I) an d eccentricity Gi pq {e) functions, respectively 
(Kaula 1966). Note that one of the major scientific 
goals of the Juno mission is a greatly improved deter- 
mination of just the harmonic coefficients of the jovian 
gravity potential; for the present-day values of the zon- 
als 41 for i = 2, 3, 4, 6 sec Table 22. According to Ander- 
son et al. (2004), it might be possible to determine the 
first three even zonals with an accuracy of 10 -9 and the 
other ones up to I = 30 at a 10 -8 level. Concerning J2, 
this would be an improvement of two orders of magni- 
tude with respect to Table 22, while the improvements 
in J4 and Jq would be of about three orders of magni- 
tude. By using the results we are going to present below 



The Jupiter gravity field is essentially determined by the Pio- 
neer 11 flyby at 1.6i?,j up (Anderson 1976); Voyager added little, 
and Galileo, which never got close to Jupiter, added nothing. 
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Table 22 Zonal harmonics of the Jupiter's gravity field according to the JUP230 orbit solution (Jacobson 2003). They 
are in unit of 10 6 . 

J2 J3 J a J& 

14696.43 ± 0.21 -0.64 ± 0.90 -587.14 ± 1.68 34.25 ±5.22 



for the long-period node and pericenter precessions, it 
can be shown that determining the low degree zonals 
at the 10 -9 level of accuracy translates into an accu- 
racy of the order of 0.5 — 1 mas yr _1 in fl and ci;, thus 
confirming the expectations of the previous Section. 
Long-period terms fulfill the condition 



£ - 2p + q = 

which results in 42 (Kaula 1966) 
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(69) 



where W, 



eop{2p-e) 



are trigonometric functions having 
the pericentre as their argument. Contrary to the small 
eccentricity satellites like the LAGEOS ones previously 
examined, in this case we will be forced to keep all 
the terms of order 0(e k ) with k > 2 in computing the 
eccentricity functions for given pairs of £ and p. More- 
over, since all the non-zero eccentricity and inclination 
functions for a given degree £ are needed, we have to 
consider all the non- vanishing terms with < p < £. 
First, we will extend our calculations to the even zon- 
als so far determined, i.e. £ = 2,4,6. In this case W 
reads 



W £0p{ 2p-i) = cos[(£ - 2p)d)t] = cos(qujt). 
It must be noted that the terms with 
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42 Note that 1(1 + 1) in eq. (69) corrects a wrong 6 in (16), pag. 
556 of Iorio (2010f). 



i.e. W, 
with 



1, yield secular precessions, while those 



q = 2p - £ ^ 



(72) 



induce harmonic signals with circular frequencies — qui. 

For degree £ = 2 only secular precessions occur and 
the non- vanishing inclination and eccentricity functions 
and their derivatives are 43 
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43 Here and in the following cl = cos/, si = sin/. Note that eq. 
(76) corrects a typo in (23), pag. 556 of Iorio (2010f). 

44 Notc that cq. (83) and eq. (84) correct typos in (30) and (31), 



p. 557 of Iorio (2010f). 
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where in addition to secular terms, also harmonic sig- 
nals with the frequencies ±2w are present. 

For £ = 6 the inclination and eccentricity functions, 
along with their derivatives read 45 
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Let us now focus on the action of the odd 
3, 5, 7, ... ) zonal (m = 0) harmonics. In this case 
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so that only harmonic terms exist for q 7^ 0. 
Finally, for I = 3 we find 46 
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Thus, we have long-period effects varying with fre- 
quencies ±w. 

We note that the long-period even and odd zonal 
harmonic terms can be approximated by secular pre- 
cessions with a high level of accuracy over the expected 
1-yr lifetime T of the Juno mission because the period 
of Junos's perijove is of the order of w 500 yr, i.e. 
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In addition to the secular rates, also harmonic signals 
with frequencies ±4w, ±2w do occur. 



sm LO 



sin(wt + w ) ~ sinw , < t < T. (103) 



Thus, the choice of the initial condition luq will be cru- 
cial in determining the impact of the long-period terms. 

Now, it would be possible, in principle, to use the 
node of Juno to measure the gravitomagnetic effect. 
Indeed, the Lense-Thirring node precession is indepen- 
dent of I, while all the zonal precessions of f2 vanish for 
I = 90 dcg. This is not the case for the perijove and the 
mean anomaly, but they are not affected by the gravit- 
omagnetic force for I = 90 deg. In reality, the situation 
will be less favorable because of the unavoidable orbit 
injection errors which will induce some departures SI 
of Juno's orbital plane from the ideal polar configura- 
tion. Thus, unwanted, corrupting node zonal secular 
precessions will appear; their mis-modeling due to the 
uncertainties S Ji may swamp the recovery of the Lense- 
Thirring effect if their determination by Juno will not 
be accurate enough. Note that there is no risk of some 
sort of a-priori "imprint" effect of the Lense-Thirring 



45 Note that eq. (94), eq. (95) and eq. (96) correct typos in (41), 
(42), (43), p. 557 of Iorio (2010f). 



46 Note that eq. (101) corrects a typo in (48), p. 557 of Iorio 
(2010f). 
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Node il only (oi = 90 deg) 




Fig. 12 Systematic percent bias on the Lense-Thirring node 
precession induced by the mis-modeling in the zonals J2 (red), 
J3 (blue), J4 (green), J$ (light blue) according to Table 22 for 89 
deg < I < 91 deg and luq = 90 deg. 



effect itself on the values of the zonals retrieved from 
the Juno's perijove motion because the gravitomagnctic 
pericenter precession vanishes for polar orbits. 

By assuming the values quoted in Table 22 for the 
uncertainties SJi, I = 2,3,4,6, let us see what is the 
impact of an imperfect polar orbital geometry on the 
nodal Lense-Thirring precessions. The results are de- 
picted in Figure 12 for each degree I separately; the ini- 
tial condition ojq = 90 deg has been used. It should be 
noted that, in view of the likely correlations among the 
determined zonals, a realistic upper bound of the total 
bias due to them can be computed by taking the linear 
sum of each mis-modeled terms. The major source of 
bias is the so far poorly known Jq; an improvement of 
four orders of magnitude, which sounds rather unlikely 
to be obtained even with Juno (Anderson et al. 2004), 
would be required to push its aliasing effect at a per- 
cent level of the Lense-Thirring effect. The situation for 
the other zonals is more favorable; J4 should be known 
better than now by a factor 1000, which is, instead, a 
realistic goal according to Anderson et al. (2004). Thus, 
we conclude that a nearly-polar orbit being 1 deg off the 
ideal 90 deg case would likely prevent a measurement 
of the gravitomagnetic node precession at a decent level 
of accuracy. 

Thanks to the high eccentricity of Juno's orbit, also 
the perijove and the mean anomaly are well defined, so 
that they can be used in a suitable way to remove the 
bias of Jq and J 2 - Let us write 



SSI = tt. 2 5J 2 + Sl. 6 8J 6 + ^n LT + Aj 



SM = M. 2 SJ 2 +M,eSJ 6 + A M , 



(104) 



(105) 



(106) 



where denotes some of Obscrvcd-minus-Calculatcd 
(O — C) quantity for the rate of the Keplerian element 'J 
which accounts for every unmodeled/mis-modeled dy- 
namical effect. It may be, for example, a correction to 
the modeled precessions to be phenomenologically esti- 
mated as a solve-for parameter of a global fit of Juno's 
data as done by Pitjeva (2005a) with the planetary peri- 
helia, or it could be a computed time-series of 47 "resid- 
uals" of ^ by suitably overlapping orbital arcs. The 
gravitomagnetic force should be purposely not modeled 
in order to be fully present in 5^ . The parameter [x is 48 
1 in GTR and in Newtonian mechanics and accounts 
for the Lense-Thirring effect. The A terms include all 
the other systematic errors like the precessions induced 
by the mis-modeled parts of the second even zonal har- 
monic SJ4 and the first odd zonal harmonic the 
mis-modeling due to the uncertainty in Jupiter's GM, 
etc. By solving for fi one obtains 



Sfl + piSus + p 2 SM = Ctur + Pi^lt + A, 



with 



Pi = 



P2 = 



M. e £1.2 - Q.6 M.2 

W.6 M.2 - M.6 UJ.2 

^.6 ^.2 — ^.6 ".2 
M.2 - M. 6 UJ.2 



W.6 



(107) 



(108) 



(109) 



Eq. (107) is designed, by construction, to single out 
the combined Lense-Thirring precessions and to cancel 
the combined secular 49 precessions due to J2 and Jq 
along with their mis-modeling. However, it is affected 
by A which acts as a systematic bias on the Lense- 
Thirring signal of interest. A globally includes the mis- 
modeled part of the combined precessions induced by 
J3 and J4; the sources of uncertainty reside in J3 and 
J4 themselves and in Jupiter's GM through the mean 
motion n which enters the coefficients Q.e , uj.e, M.e. 

In Figure 13 the impact of the mis-modeling in J3 
and J 4 for uj = 90 deg is depicted. It corrects Fig. 2, 
p. 558 of Iorio (2010f); incidentally, Figure 13 yields 
smaller percent uncertainties. In Figure 14, which 
amends Fig. 3, p. 559 of Iorio (2010f), we use ujq = 
deg. In this case the situation is much more favorable 
because for a total departure of ±1 deg from I = 90 
deg, an improvement of only two orders of magnitude 
in J3, which is, today, still compatible with zero, and J4 

47 Since the Keplerian elements are not directly measurable quan- 
tities, we use here the term "residual" in an improper sense. 

48 It is not one of the standard PPN parameters, but it can be 
expressed in terms of 7 as fj, = (1 + 7)/2. 

49 We include in them also the long-period harmonic terms for 
the reasons explained before. 
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li + Ciw + c 2 M (ajQ = 90 deg) 




Fig. 13 Systematic percent bias on the Lense-Thirring pre- 
cessions, combined according to eq. (107), induced by the mis- 
modeling in the uncanceled zonals J3 (dash-dotted red line), J4 
(dotted blue line) according to Table 22 for 89 deg < / < 91 deg 
and ujq = 90 deg. 

H + CtO) + c 2 M (oi = deg) 
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Fig. 14 Systematic percent bias on the Lense-Thirring pre- 
cessions, combined according to eq. (107), induced by the mis- 
modeling in the uncanceled zonals J3 (dash-dotted red line), J4 
(dotted blue line) according to Table 22 for 89 deg < / < 91 deg 
and uq = deg. 



would be needed to reach the percent level; let us recall 
that the expected improvement in J 4 with respect to 
the results by Jacobson (2003) is about three orders of 
magnitude (Anderson et al. 2004). Note that a value 
of u)q far from 90 deg is preferable to minimize the per- 
turbations. 

Another potential source of systematic errors is 
Jupiter's GM whose uncertainty 5(GM) indirectly af- 
fects eq. (107) through the Keplerian mean motion 
n entering the uncanceled J3 and J4 combined pre- 
cessions; Sn is also present via the mean anomaly it- 
self. However, it turns out that it is of no concern 
because, according to the present-day level of relative 
uncertainty (Jacobson 2003) 



5{GM) 
GM 



= 1.6 x 10~ 8 , 



(110) 



its impact on the combined Lense-Thirring precessions 
is well below the percent level. 

6.2.2 A numerical approach 

Also in this case we followed an alternative approach 
based on preliminary numerical simulations. We inves- 
tigated the impact of the uncertainties in the first two 
jovian even zonals on a Juno's single six-hours pass by 
numerically simulating the probe's Doppler range-rate 
signals due to 8J2 and 6J4. By assuming for them val- 
ues as large as 2 x 10~ 10 and 3 x 10 -10 , respectively, it 
turns out that the maximum Doppler shifts arc roughly 
1 — 1.5 jjLva s . Moreover, and more importantly, the 
time-dependent patterns of the even zonals' Doppler 
signals are quite different from the Lense-Thirring one 
removing the risk of an insidious mimicking bias. An- 
other encouraging fact is that such simulations indicate 
that an inclination of even 91 deg would not compro- 
mise the recovery of the gravitomagnctic signal of in- 
terest. 

6.3 Summary 

We first explored the possibility of a high accuracy 
measurement of the Lense-Thirring effect by perform- 
ing analytical calculations and interpreting them in a 
rather conservative fashion. A meter-level accuracy in 
determining the jovicentric orbit of the Juno space- 
craft should not be an unrealistic goal to be reached. 
Equivalcntly, the gravitomagnctic node precession of 
Juno amounts to 68.5 mas yr , while the accuracy 
in measuring its node and pcrijovc precessions should 
be of the order of 0.5 — 1 mas yr -1 , given the ex- 
pected improvements in our knowledge of the depar- 
ture of the jovian gravitational field from spherical sym- 
metry. If Juno's orbit were perfectly polar, the long- 
period node precessions induced by the zonal harmonics 
Ji, £ = 2, 3, 4, 6, . . . of the non-spherical jovian gravita- 
tional potential would vanish, thus removing a major 
source of systematic alias on the Lense-Thirring secular 
precession. In reality, unavoidable orbit injection errors 
will displace the orbital plane of Juno from the ideal po- 
lar geometry; as a consequence, the mis-modeled part 
of the node zonal precessions would overwhelm the rel- 
ativistic signal for just SI = ±1 deg, in spite of the 
expected improvements in Ji,£ = 2,4,6 of three or- 
ders of magnitude. A suitable linear combination of 
the node tt, the perijove w and the mean anomaly M. 
will allow to cancel out the effects of J2 and Jq; the 
remaining uncanceled J3 and J4 will have an impact on 
the Lense-Thirring combined precessions which should 
be reduced down to the percent level or better by the 
improved low-degree zonals. 
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Instead of looking at the cumulative, secular effects 
over the entire duration of the mission, we also followed 
an alternative approach by looking at single Doppler 
range-rate measurements over time spans six hours 
long centered on the the probe's closest approaches to 
Jupiter; it turned out that, in this way, the perspec- 
tives are even more favorable. We numerically sim- 
ulated the characteristic Lense-Thirring pattern for a 
single science pass by finding a maximum value of the 
order of hundreds /mi s _1 , while the expected precision 
level in Juno's Doppler measurements is of the order 
of one /zm s _1 . Thus, by exploiting about 25 of the 
planned 33 total passes of the mission it would be possi- 
ble to reach a measurement accuracy below the percent 
level. We repeated our numerical analysis also for the 
Doppler range-rate signals of Ji and J4 by finding quite 
different patterns with respect to the gravitomagnctic 
one; moreover, for a level of mismodcling of the order 
of 2 — 3 x 10~ 10 in such zonals the maximum value 
of their biasing Doppler signals is about 1 — 1.5 /mi 
s _1 . Our numerical analysis also shows that a depar- 
ture from the nominal polar orbital geometry as large 
as 1 deg would not compromise the successful outcome 
of the measurement of interest, contrary to the con- 
servative conclusions of our analytical analysis. Thus, 
this approach suggests that there is not a high corre- 
lation between the Lense-Thirring parameter and the 
jovian gravity held parameters, although a covariance 
analysis would be needed to prove it. However, such a 
covariance analysis is outside the scope of the present 
paper. 

In conclusion, the potential error in the proposed 
Juno Lense-Thirring measurement is between 0.2 and 
5 percent. Conversely, if one assumes the existence 
of gravitomagnctism as predicted by general relativity, 
the proposed measurement can also be considered as 
a direct, dynamical determination of the jovian proper 
angular momentum S by means of the Lense-Thirring 
effect at the percent level. 

7 A note on gravitational waves: a new 
window into the Universe 

In the Introduction we emphasized that, as gravitomag- 
netic effects are usually discussed in the linearized ap- 
proach of GTR, they also have a strong connection with 
GWs within various theories of gravity (Iorio and Corda 
2009, 2010; Corda et al. 2010). 

In fact, interferometric GW detectors are operative 
at the present time (Giazotto 2008; Iorio and Corda 
2010). A direct detection will be a historic confirma- 
tion of the indirect evidence of the existence of GWs 



by Hulse and Taylor (1975). The realization of a GW 
astronomy, by giving a significant amount of new infor- 
mation, will be a cornerstone for a better understanding 
of gravitational physics too. Detectors for GWs will 
be important for a better knowledge of the Universe 
(Giazotto 2008) and also because the interferometric 
GW detection will be the definitive test for GTR or, 
alternatively, a strong endorsement for Extended The- 
ories of Gravity (Corda 2009). In fact, if advanced 
projects on the detection of GWs improve their sensitiv- 
ity, allowing the scientific community to perform a GW 
astronomy, accurate angle- and frequency-dependent 
response functions of interferometers for GWs arising 
from various theories of gravity will permit to discrimi- 
nate among GTR and extended theories of gravity like 
Scalar Tensor Gravity, String Theory and f(R) The- 
ories. This ultimate test will work because standard 
GTR admits only two polarizations for GWs, while in 
all extended theories the number of polarization states 
is, at least, three; see Corda (2009) for details. On 
the other hand, the discovery of GW emission by the 
compact binary system composed of two neutron stars 
PSR1913+16 (Hulse and Taylor 1975) has been, for 
both of theoretician and experimental physicists work- 
ing in this field, the fundamental thrust allowing to 
reach the extremely sophisticated technology needed 
for investigating this field of research (Iorio and Corda 
2010). 

GWs arise from Einsteinian GTR (Einstein 1915) 
and are weak perturbations of the spacetime curvature 
which travel at light speed (Einstein 1916, 1918). As 
GWs arc quadrupolar waves, only asymmetric astro- 
physical sources can emit them (Iorio and Corda 2010). 
The most efficient ones are very dense and massive co- 
alescing binary systems, like neutron stars and black 
holes. Instead, a single rotating pulsar, even if dense 
and massive, can only rely on spherical asymmetries, 
usually very small (Iorio and Corda 2010). Supernova 
explosions could have, in principle, relevant asymme- 
tries, being potential sources of gravitational waves (Gi- 
azotto 2008; Iorio and Corda 2010). The most im- 
portant cosmological source of GWs is the so called 
stochastic background of relic GWs which, together 
with the Cosmic Background Radiation (CBR), carries 
a huge amount of information on the early stages of 
the Universe evolution (Smoot and Steinhardt 1993). 
This reason is a rather hopeful one (Smoot and Stein- 
hardt 1993; Allen 1997; Corda 2010). As gravitation 
is the weakest of the four known fundamental inter- 
actions, the small-scale perturbations of the gravita- 
tional field decoupled from the evolution of the rest 
of the universe at very early times. Currently, the 
most detailed view of the early universe comes from 
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the microwave background radiation, which decoupled 
from matter about 300.000 years after the big bang, 
and gives us an accurate picture of the universe at 
this very early time (Smoot and Steinhardt 1993; Allen 
1997; Corda 2010). Some simple computations by Allen 
(1997) showed that if a background of relic GWs is ob- 
served, they will carry a picture of the universe as it 
was about 10~ 22 s after the initial singularity. This 
would represent a tremendous step forward in our un- 
derstanding of the primordial universe. The existence 
of a relic stochastic background of GWs is a conse- 
quence of very general assumptions which arise from 
an interplay between basic principles of classical theo- 
ries of gravity and of quantum field theory (Smoot and 
Steinhardt 1993; Allen 1997; Corda 2010). The strong 
variations of the gravitational field in the primordial 
universe amplified the zero-point quantum oscillations 
and produced relic GWs. This model derives from the 
inflationary framework of the early universe (Lyth and 
Liddle 2009), which is fine tuned with the WMAP data 
on the CBR. In fact, an exponential inflation with a 
spectral index of f=a 1 agrees with inflationary theories 
(Spergel et al. 2003). The inflationary scenario predicts 
cosmological models in which the Universe undergoes a 
brief phase of a very rapid expansion in early times 
(Lyth and Liddle 2009). Such an expansion could be 
either a power-law or exponential in time (Lyth and 
Liddle 2009) and provides solutions to the horizon and 
flatness problems. Inflation also contains a mechanism 
which creates perturbations in all fields (Allen 1997; 
Corda 2010; Lyth and Liddle 2009). This mechanism 
also yields a distinctive spectrum of relic GWs (Allen 
1997; Corda 2010). The GW perturbations arise from 
the uncertainty principle and the spectrum of relic GWs 
is generated from the adiabatically-amplificd zero-point 
fluctuations (Allen 1997; Corda 2010). 

Regarding the potential GW detection, recalling 
some historical notes is in order. In 1957 F.A.E. Pirani, 
a member of Hermann Bondi's research group, pro- 
posed the geodesic deviation equation as a theoretical 
foundation for designing a practical GW detector (Pi- 
rani 1957). In 1959, J. Weber studied a detector that, 
in principle, might be able to measure displacements 
smaller than the size of the nucleus (Weber 1959). He 
realized an experiment using a large suspended bar of 
aluminum. Such a bar had a high resonant Q at a fre- 
quency of about 1 kHz. After this, in 1960, he tried to 
test the general relativistic prediction of gravitational 
waves from strong gravity collisions (Weber 1960). We- 
ber performed further analyses in 1969, by claiming ev- 
idence for the observation of gravitational waves from 
two bars separated by 1000 km (Weber 1969). Such 
evidence was based on coincident signals which Weber 



claimed to be detected by the two bars. In order to 
confirm these observations, he also proposed to detect 
gravitational waves by using laser interferometers (We- 
ber 1969). In fact, all the modern detectors can be 
considered to originate from early Weber's studies (Gi- 
azotto 2008; Iorio and Corda 2010). 

Currently, five cryogenic bar detectors have been 
built to work at very low temperatures (< 4 K): Ex- 
plorer at CERN, Nautilus at Frascati INFN National 
Laboratory, Auriga at Legnano National Laboratory, 
Allegro at Louisiana State University and Niobe in 
Perth (Giazotto 2008; Iorio and Corda 2010). There are 
also two spherical detectors, i.e. the Mario Schenberg, 
which operates in San Paolo (Brazil) and the Mini- 
GRAIL, which is trying to detect GWs at the Kamcr- 
lingh Onncs Laboratory of Leiden University; see Gi- 
azotto (2008) and Iorio and Corda (2010). Spherical 
detectors are quite important for the potential detec- 
tion of the scalar component of GWs that is admitted 
by Extended Theories of Gravity (Corda 2007a). In 
the case of interferometric detectors, two mirrors are 
used like free falling masses. Such mirrors are sepa- 
rated by by 3 km for Virgo and 4 km for LIGO. Thus, 
the GW tidal forces in interferometers are expected to 
be several order of magnitude larger than in bar detec- 
tors. Differently from bars, interferometers like LIGO 
have very large bandwidth (10-10000 Hz) because mir- 
rors are suspended to pendulums having resonance in 
the Hz region. In this way, above such a resonance fre- 
quency, mirrors works, in a good approximation, like 
freely falling masses in the horizontal plane (Giazotto 
2008; Iorio and Corda 2010). 

Recently, starting from the analysis in Baskaran and 
Grishchuk (2004), some papers have shown the impor- 
tance of the gravitomagnetic effects in the framework of 
the GW detection, too; see Iorio and Corda (2010) and 
Corda (2007b) and references therein. In fact, the so- 
called magnetic components of GWs have to be taken 
into account in the context of the total response func- 
tions of interferometers for GWs propagating from arbi- 
trary directions (Baskaran and Grishchuk 2004; Corda 
2007b; Iorio and Corda 2010). An extended analysis, 
which has been carefully reviewed in the recent work 
(Iorio and Corda 2010), showed that such a magnetic 
component becomes particularly important in the high- 
frequency portion of the range of ground based inter- 
ferometers for GWs which arises from standard GTR. 
The magnetic component has been extended also to 
GWs arising from scalar-tensor gravity (Corda et al. 
2010; Iorio and Corda 2010). Such studies showed that 
by considering only the low-frequency approximation 
of the electric contribution and thereby neglecting the 
magnetic contribution, a portion of about 15% of the 
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signal could be, in principle, lost in the case of Scalar 
Tensor Gravity too (Corda et al. 2010; Iorio and Corda 
2010), in close analogy with the standard case of GTR 
(Baskaran and Grishchuk 2004; Corda 2010). 

It is important to give a physical/mathematical ex- 
planation as to why the magnetic component of the 
gravitational field of a GW becomes important at high 
frequency. As interferometric GWs detection is per- 
formed in a laboratory environment on Earth, the co- 
ordinate system in which the spacetime is locally flat 
is typically used (Iorio and Corda 2010) and the dis- 
tance between any two points is given simply by the 
difference in their coordinates in the sense of Newto- 
nian physics. In this frame, called the frame of the 
local observer, GWs manifest themselves by exerting 
tidal forces on the masses (the mirror and the beam- 
splitter in the case of an interferometer) . In the follow- 
ing, we work with G = 1, c = 1 and h = 1 and we call 
h+(ttt+ztt) and h x (Ut+ztt) the weak perturbations due 
to the + and the x polarizations which are expressed 
in terms of synchronous coordinates t t t,xtt,ytt, z tt in 
the transverse-traceless (TT) gauge (Iorio and Corda 
2010). In this way, the most general GW propagating 
in the Ztt direction can be written in terms of plane 
monochromatic waves (Iorio and Corda 2010; Baskaran 
and Grishchuk 2004) 



zu) = h+(tu + Ztt)e$ '+ 
+ h x (ttt + tkt)e$ = 
= h + oexpioj(t u + ztt)e$/ + 



(111) 



(x) 



+ h x0 exp iu>(t tt + Zu)eii V 
with the corresponding line element 

ds 2 = dt 2 t - dz 2 tt - (1 + h + )dx% - (1 - h + )dyf t - 



- 2h x dx tt dx t t- 

(112) 

The wordlines Xtt, ytt, z u = const, are timelike geodesies 
representing the histories of free test masses (Iorio 
and Corda 2010). The coordinate transformation 
x a = x a (xt t ) from the TT coordinates to the frame 
of the local observer is given by (Iorio and Corda 2010; 



Baskaran and Grishchuk 2004) 
t = t tt + \{x 2 tt - y 2 t )h + - \x tt ytth x 



Xtt 



+ \x tt h + - \y u h x + \x u z t th-\ 



hyttztthx 



+ \y u h+ - hxuhy + \y u z t th+- 



(113) 



\xttZtthy, 

Ztt - z( x tt - Vtt)h+ + \xttytth, 



and 

/' + 



dh + 

dt ■ 

dh x 
dt ■ 



(114) 



The coefficients of this transformation (components of 
the metric and its first time derivative) are taken along 
the central wordline of the local observer (Iorio and 
Corda 2010; Baskaran and Grishchuk 2004). It is well 
known that the linear and quadratic terms, as pow- 
ers of unambiguously determined by the con- 
ditions of the frame of the local observer, while the 
cubic and higher-order corrections are not determined 
by these conditions (Iorio and Corda 2010; Baskaran 
and Grishchuk 2004). Thus, at high-frequencies, the 
expansion in terms of higher-order corrections breaks 
down (Iorio and Corda 2010; Baskaran and Grishchuk 
2004). Considering a free mass riding on a timelike 
geodesic (x = h, y = I2, z = £3 ) (Iorio and Corda 
2010; Baskaran and Grishchuk 2004), eq. (113) defines 
the motion of this mass with respect to the introduced 
frame of the local observer. In concrete terms one gets 

x(t) = h + \[hh+{t)-hhx{t)] + \hhh+{t)+ 

+ \hhhx{t) 

y(t) = h-\[kK{t) + hhx{t)]-\hhh+{t)+ 

+ \hhhx{t) 
z(t) = h-\{l 2 -l 2 )h + {t) + 2hl 2 h y {t). 

(115) 

In the absence of GWs the position of the mass is 
(^1^2,^3) and the effect of the GW is to induce oscil- 
lations of the mass. Thus, in general, from eq. (115) 
all three components of motion are present (Iorio and 
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Corda 2010; Baskaran and Grishchuk 2004). Neglect- 
ing the terms with h + and h x in eq. (115), the "tra- 
ditional" equations for the mass motion are obtained 
(Iorio and Corda 2010; Baskaran and Grishchuk 2004) 

x(t) = h + \[hh + (t)-l 2 h x (t)] 

y(t) = l2-\[hh + {t)+hh x {tj\ (116) 

z{t) = h. 

Clearly, this is the analogon of the electric component 
of motion in electrodynamics (Iorio and Corda 2010; 
Baskaran and Grishchuk 2004), while the equations 

x{t) = h + \hhh + {t) + \hhh x (t) 

y(t) = h-\hhh+{t) + \hhh x (t) (117) 

*(*) = h-\{l 2 1 -ll)h + {t) + 2l 1 l 2 h x (t), 

are the analogues of the magnetic component of motion. 
One could think that the presence of these "magnetic" 
components is a "frame artefact" due to the transfor- 
mation of eq. (113), but in Section 4 of Baskaran and 
Grishchuk (2004) eq. (115) has been directly obtained 
from the geodesic deviation equation too, thus the mag- 
netic components have a real physical significance. The 
fundamental point of Iorio and Corda (2010); Baskaran 
and Grishchuk (2004) is that the "magnetic" compo- 
nents become important when the frequency of the 
wave increases but only in the low- frequency regime. 
This can be understood directly from eq. (115). In 
fact, by using eq. (Ill) and eq. (113), it turns out that 
eq. (115) becomes 

x(t) = l 1 + ±[l 1 h + (t)-l 2 h x {t)}+ 

+ \hhuh+{t-^) + \hhuh x (t-l) 
y{t) = l 2 -\[hh + {t)+hh x {t)]- 

(118) 

- \hhuh+{t - § ) + \hhujh x [t - f ) 

*(*) = i3-^(i?-£)wM*-5) + 

+ 2hl 2 uh x (t-%). 

Thus, the terms with h+ and h x in eq. (115) can be ne- 
glected only when the wavelength goes to infinity, while, 
at high-frequencies, the expansion in terms of colilj cor- 
rections, with i, j = 1,2,3, breaks down (Iorio and 
Corda 2010; Baskaran and Grishchuk 2004). The at- 
tentive reader could be surprised that gravitomagnetic 



effects in the field of a GW have been ignored in the 
past. Actually, the key point is that realizing the first 
direct detection has been always the most important 
goal for scientists in this research field. To obtain such a 
first detection, the low frequency approximation is con- 
sidered sufficient, i.e. only the "electric" component 
is needed for ground based interferometers (Baskaran 
and Grishchuk 2004). On the other hand, if the scien- 
tific community aimes to realize a concrete GW astron- 
omy, the above discussion shows that gravitomagnetic 
effects in the field of a GW have to be taken in due 
account. Finally, we recall that GM effects arc im- 
portant in the process of producing GWs too. Indeed, 
Capozzicllo ct al. (2010) studied corrections to the rela- 
tivistic orbits by considering high order approximations 
induced by GM effects for very massive binary systems 
(e.g. a very massive black hole on which a stellar object 
is spiralling in). Such corrections have been explained 
by taking into account "magnetic" components in the 
weak field limit of the gravitational field. New nuta- 
tion effects of order c~ 3 have been found. These effects 
work beside the standard periastron corrections and af- 
fect the GW emission and the gravitational waveforms 
due by the orbital motion of the two massive bodies 
(Capozziello et al. 2010). The GM corrections emerge 
as soon as matter-current densities and vector gravi- 
tational potentials cannot be discarded into dynamics 
(Capozziello et al. 2010). GWs emitted through mas- 
sive binary systems have been studied in the quadrupole 
approximation (Capozziello et al. 2010). 

8 Conclusions 

Measuring gravitomagnctism, and, in particular, the 
Lense-Thirring effect is a challenging and difficult en- 
terprize. In this review we have shown how the de- 
tection of such a general relativistic signal is closely 
related to a number of other competing dynamical fea- 
tures of motion, so that a genuine interdisciplinary at- 
titude is required. Looking for the prediction by Lense 
and Thirring may help in shedding light on different 
aspects of the solar system scenarios in which it is at- 
tempted to be measured. Conversely, a better knowl- 
edge of the forces acting on dynamical systems may be 
of great help in reliably and accurately measuring the 
gravitomagnetic field itself. 

The current tests with the LAGEOS satellites in 
the gravitational field of the Earth, which represent 
the first implemented attempts to measure the Lense- 
Thirring effect, should be repeated by directly mod- 
eling the gravitomagnetic force and explicitly solving 
for it in the data processing. Moreover, they should 
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be complemented by a re-analysis of the data sets of 
the currently ongoing spacecraft-based missions dedi- 
cated to accurately measuring the global gravitational 
field of the Earth like GRACE taking into account gen- 
eral relativity itself; new global solutions in which one 
or more parameters accounting for the Lense-Thirring 
effect are solved-for should be produced. For the mo- 
ment, a conservative evaluation of the systematic bias 
in the LAGEOS-LAGEOS II tests due to the imperfect 
knowledge of the classical part of the Earth's gravita- 
tional field points toward a 20 — 30% level; more op- 
timistic evaluations yield a 10 — 15% uncertainty. It 
is unclear if the future LARES mission will be able to 
measure the Lense-Thirring precessions at the repeat- 
edly claimed « 1% accuracy. Indeed, its orbital con- 
figuration, different from that originally proposed, may 
enhance some competing dynamical effects of gravita- 
tional and non-gravitational origin at a level whose un- 
certainty is difficult to be realistically evaluated. 

The present-day level of accuracy in the orbit deter- 
mination of the inner planets of the solar system has 
basically reached the magnitude of the Lense-Thirring 
precessions. The most important systematic error due 
to mis-modeled dynamical effects is associtated with 
the first even zonal harmonic of the non-spherical grav- 
itational field of the Sun, which particularly affects Mer- 
cury. Anyway, an accurate measurement of the Sun's 
oblateness is one of the goals of future spacecraft-based 
missions. The analysis of more accurate ranging data 
to present and future spacecraft like, e.g., BepiColombo 
and, especially, future interplanetary laser ranging de- 
vices should drastically improve the situation allow- 
ing for a reasonably accurate measurement of the solar 
gravitomagnctic field. 

After the preliminary attempts with the Mars Global 
Surveyor probe which have paved the way, also Mars 
may become suitable for reliably measuring the Lense- 
Thirring orbital precessions in the near future. To this 
aim, the forthcoming Phobos-Grunt mission, scheduled 
for launch in late 2011 or early 2012, will improve our 
knowledge of some key physical properties of Mars and 
of the orbit of its satellite Phobos. 

The Juno mission to Jupiter, scheduled for launch in 
2011, may yield the opportunity for an accurate mea- 
surement of the Jovian angular momentum through the 
Lense-Thirring effect itself. 

Finally, we emphasized the relevant connection be- 
tween gravitomagnetic effects and GWs. In fact, the 
scientific community expects the first direct detection of 
GWs within the next years. The importance of gravit- 
omagnctic effects in the framework of GWs increases in 
the high-frequency portion of the range of ground based 
interferometers for both of GWs which arise from stan- 



dard GTR as well as those coming from scalar-tensor 
gravity. If one neglects the magnetic contribution con- 
sidering only the low-frequency approximation of the 
electric contribution, a portion of about 15% of the 
signal could be lost by interferometric detectors. Thus, 
a carefully analysis of the "magnetic" -type contribution 
of GWs is needed in order to establish a GW astronomy 
within the next years. 
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